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1. Introduction 

Geometric scattering theory, as espoused in jl9| , is the study of natural operators 
such as the Laplacian on non-compact Ricmannian manifolds X with controlled 
asymptotic geometry using geometrically-informed, fully (i.e. to the extent it is 
possible) microlocal methods. The goals of this subject include the construction of 
an analytically useful compactification of X and the definition of an appropriate 
class of pseudodifferential operators containing the resolvent of the Laplacian, the 
Schwartz kernel of which is a particularly simple distribution on this compactifica- 
tion. In this paper we examine the case where X is a product of asymptotically 
hyperbolic (or conformally compact, as they are often called) spaces from this point 
of view. This is intended both as an initial application of these methods to higher 
rank symmetric spaces and their geometric generalizations, and also as a relatively 
simple (although still surprisingly complicated) example which should provide a 
guide for what to expect in further development in this area. 

The results here tend to be notationally complicated, so for the purposes of this 
introduction we state a simple, yet representative, result concerning the asymptotics 
of the resolvent applied to a Schwartz function. Suppose that (Mj,gj), j = 1, 2, are 
(conformally compact) asymptotically hyperbolic, dimMj = kj + 1, Hj = ~A gj . 
positive, 4>ji are the L 2 eigenf unctions of Hj with eigenvalue Xji < fcJ/4. The usual 
compactification of Mj as a C°° manifold with boundary is denoted Mj, and has 
boundary defining function Xj. Adjoining pj — — 1/logXj to the smooth structure 
of Mj yields the logarithmic blow up (Mj)i og . Let H be the sum of the Laplacians 
from the two factors: H = Hi <E> Id + Id^i^, and R(p) the resolvent, (H — 1 . 
Define 

(1.1) X = [(Mi)io g x (M 2 ) log ;a(M 1 ) log x 9(M 2 ) log ], 

which is a resolution of X = M\ x M%. By definition, each pj = — 1/logXj is 



smooth on (Mj)i og and thus p, where p 1 = y p 1 2 + p 2 2 is smooth on X 
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Theorem. (See Theorem |7^| .j Suppose f e C°°(X), R9(i> fc 2 /4 = (fc? + fel)/ 4 - 
Then on X: 



(1.2) 



- iO)/ = 4 l/2 X2 2/2 exp(-i v /^~fc2/4/ j o).g 
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g polyhomogeneous on X, gu and gu are polyhomogeneous on Mi and Mi respec- 
tively. 

The leading term of each part in the asymptotics can be described explicitly. 
Since 0™ decays like x] 3 t 2+ v^ ; the first term dominates the other two in 



the interior of the front face of the blow-up (1.1); the second term, involving the 
eigenfunctions of Hi, is only comparable to it at the lift of Mi x dM 2 , while the 
third term, involving the eigenfunctions of Hi, is only comparable to it at the lift 
of dMi x M 2 . 

Similar results are valid outside the continuous spectrum, where the exponents 
are not pure imaginary, but care needs to be taken as to which eigenvalue terms 
appear in the asymptotics. Moreover, additional singularities appear along the 
front face of X when fi is real, below the spectrum of H , which, however, can be 
resolved by further re al b low-ups (or understood as a type of Legendre singularity). 



We refer to Theorem |6.3| for the detailed statement of the results in that case. We 
deduce similar results for the resolvent kernel, which we use in turn to analyze 
the Martin compactification Xmo{ Mi x Mi. While the latter behaves (nearly) 
as expected when Hi and Hi have no L 2 eigenvalues, it experiences a substantial 
collapse in the presence of such eigenvalues. 

Theorem. There is a natural continuous surjection X — > Xm- Moreover, the 
following hold. 

(i) If neither Hi nor Hi have L 2 eigenvalues, the rest riction of this map to a 



neighborhood of the front face of the blow-up in ( hi ) is injective. In general. 



its injectivity on the 'side faces ' of X depends on properties of the spherical 
functions, i.e. on -^\ T =oRj{k 2 + t 2 ). 
(ii) If both H\ and H 2 have L 2 eigenvalues, the Martin boundary dX m is of the 
form dM 1 U dMi U <9Mi x dM 1 x I , I an open interval. This set is naturally 
identified with a collapsed version of X, and the map dX — > dXu factors 
through it. 

There is an extensive literature on scattering on conformally compact manifolds, 
including especially the important special case of convex cocompact (and geomet- 
rically finite) hyperbolic manifolds. This contains geometric constructions of the 
resolvent, scattering matrix and generalized eigenfunctions, as well as trace formulae 
and asymptotics for the counting function for resonances. For brevity, we list only 
(H), |1, @, Hi, @, 13 for references. A parallel theory for complex hyperbolic 
manifolds and their perturbations is initiated in Q , and while not appearing explic- 
itly, this theory extends to manifolds with the asymptotic structure of quaternionic 
hyperbolic spaces and the Cayley plane. 

Concerning higher rank noncompact symmetric spaces (and in less generality, 
their quotients too), the compactification theory is now well-understood Q, as well 
as at least some aspects of the analysis of the Laplacian. However, the current 
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methods here rely heavily on the special algebraic structure of these spaces, so 
extensions of these results to even relatively modest geometric perturbations of 
these spaces are basically not understood at all. 

Alongside this is the study of the asymptotically Euclidean scattering metrics 
initiated in ]l8| , pp|] . This theory extends the extensive classical literature, but has 
led to a new and detailed understanding of quantum N-body scattering where the 
geometry has much in common with that of flats in non-compact symmetric spaces, 
see e.g. ||, f§, @ and j2jj. 

These different geometric settings have led to the understanding of diverse ana- 
lytic phenomena which may sometimes be traced to the effects of the asymptotically 
flat or negative curvature. One of the attractions of studying higher rank symmet- 
ric spaces from the point of view of geometric scattering theory is to isolate the 
specific ways in which the flat and negatively curved directions interact with one 
another and affect the analysis. The simplest setting where these sorts of effects 
might be seen is on the product X = Mi x M2, where both factors (Mj,gj) are 
conformally compact metrics. Recall that (M, g) is conformally compact if M is 
a smooth compact manifold with boundary, such that for some defining function 
x for dM, g takes the form (dx 2 + h)/x 2 , where ft, is a nonnegative smooth sym- 
metric 2-tensor which restricts to a nondegenerate metric on dM. (Note, however, 
that only the conformal class of h on dM is well-defined from g.) The prototype 
of a conformally compact manifold is hyperbolic space, and so the prototypes for 
the spaces we consider here are products of hyperbolic spaces. The reader should 
be aware, though, that Mi x M% is asymptotically like the product of hyperbolic 
spaces H™ 1 x H™ 2 (rii — dimM^) only near dM\ x dM2, but elsewhere may be a 
rather severe metric and topological perturbation of this model space. 

The basic questions we ask here concerning the product space X with metric 
g = gi + 32 concern the resolvent RxW = (A g — A) -1 of its Laplacian A g . We 
sometimes write i?x(A) as R g (X) or simply R(X). This is a holomorphic family of 
elliptic pseudodifferential operators of order —2 in the resolvent set C \ spec(A 3 ), 
which contains at least the complement of the positive real axis. The precise ex- 
tent of the spectrum is straightforward to determine from the spectra on either 
factor, but the first substantial question is to understand the behaviour of R(X) 
as A approaches the (continuous part of the) spectrum. Existence of a limit, in 
an appropriate sense, is known as the limiting absorption principle. Even better is 
the existence of a meromorphic continuation of -R(A) beyond the spectrum. Usu- 
ally, when it exists, this continuation lives on some Riemann surface covering the 
complex plane and ramified at the thresholds of the spectrum. Meromorphic con- 
tinuations of this type are known to exist for the Laplacian of conformally compact 
manifolds p"5| . Further important questions concern the geometric space which is 
a resolution of X x X obtained by a process of real blow-up, and on which the 
Schwartz kernel of the resolvent lives as a polyhomogeneous (or more generally, a 
Legendrian) distribution. Detailed knowledge of this space is central in understand- 
ing the finer properties of the resolvent, and conversely is key in its initial geometric 
construction. Geometric and analytic compactifications of the space X itself are 
quite relevant to this. There are two well-known compactifications which are closely 
related to our methods: the geodesic (also known as the conic) compactification, 
and the Martin compactification, which is defined using function theory, specifically 
the space of positive solutions of (A g — X)u = for A real and below the bottom of 
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spec(A g ). This Martin compactification is known in many instances, including for 
the product of hyperbolic spaces (|] and more recently for general symmetric spaces 
of noncompact type jjj. The 'resolvent compactification' of X x X we construct 
below, and on which the Schwartz kernel of the resolvent lives, has a structure on 
some of its hypersurface boundary faces which is inherited from the geodesic and 
Martin compactifications. 

These questions together have led to the somewhat modest goals of the present 
paper. In the next section we present a contour integral formula for Rx(X) in terms 
of the resolvents of the two factors RmjW- A related expression, written as an 
integral over the spectral measure, derived in the context in Euclidean scattering, 
appears in work of Ben-Artzi and Devinatz J3|. The representation formula here 
holds in great generality. In §3 we specialize and review the detailed structure 
of the resolvent Rm when M is a conformally compact manifold; some auxiliary 
estimates required later are also derived here. An immediate consequence of the 
representation formula of §2 is the existence of an appropriate meromorphic contin- 
uation for the resolvent Rx (A), and we describe this in §4. After this, §5 contains a 
discussion of general compactification theory, specialized to this context. The main 
work is done in §6 and §7, where we describe the asymptotics of Rx{fJ>)f when / is 
Schwartz, first when /i is in the resolvent set and then when /i is in the main sheet 
of continuous spectrum of H. The first application of this is given in §8, where 
we construct the 'resolvent double space', a resolution of X x X which carries the 
Schwartz kernel of Rx{^) in as simple a fashion as possible. Finally, in §9, we use 
the resolvent asymptotics to determine the Martin compactification of X. 

There are several features of this work to which we wish to draw particular 
attention. First, the main tool in deriving the resolvent asymptotics is stationary 
phase, which is in many respects local. The previous identification of the Martin 
compactification (at least for the product of hyperbolic spaces) in J(| relies heavily 
on global heat kernel bounds, which we feel are intrinsically more complicated. Note 
also that by using stationary phase we are taking advantage of the oscillatory nature 
of the resolvent, even when studying it for certain values of /i where it has been 
more traditional to rely on 'positivity methods' such as the maximum principle, 
the Harnack inequality, etc. Another interesting feature here is the surprisingly 
complicated way the existence of bound states, i.e. L 2 eigenvalues, for the Laplacian 
on either factor affects the asymptotics and the structure of the Martin boundary. 
Finally, we have given a detailed description of the smooth structure of the various 
compactifications we construct; this aspect is usually neglected in other discussions 
of compactification theory, but as we show, plays a significant role. 

Our intention is that the investigations here will form the basis for a more thor- 
ough investigation of the geometric scattering theory for higher rank spaces. 

The authors wish to thank Richard Melrose for helpful advice, and also Lizhen 
Ji for encouraging us to study the Martin compactification. 

2. Resolvent formula 

Let Hi, H2 be sclf-adjoint operators on Hubert spaces Vi and V2 which are 
bounded below; the precise structure of their spectra will be unimportant for the 
present. We denote inf spec(ifj) = Ao(-Hj). Now let 



H = Hi (g> Id + Id (g>H 2 
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be the self-adjoint operator on the completed tensor product space Vi®V2- Then 
inf spec(if) = X (H) = A (#i) + X Q (H 2 ). Let 

Rj(Hj) = {Hj - fij)- 1 , ; = 1,2, and - (H - m)" 1 

be the resolvents of the Hj and H, respectively. This setting has been inves- 
tigated by Ben-Artzi and Devinatz in p), mostly from the point of view of the 
limiting absorption princple, i.e. the existence of the boundary values R(p ± iO), 
M € [Ao(-ffi) + Ao(i?2),oo) on suitable weighted spaces, under the assumption that 
the limiting absorption principle holds for the Rj individually. 

One of our first goals is to show that if both Rj admit meromorphic continuations, 
then R(fj.) does as well. Later we wish to obtain precise asymptotics for the Schwartz 
kernel of R(fJ.) for \x in the resolvent set and in the continuous spectrum. In this 
section we derive a representation of R(fi) as a contour integral which will be useful 
for both of these purposes. 

So, fix m in the resolvent set C\ [Ao(i?i)+Ao(-ff2), oo), and let 7 be a parametrized 
curve in the complex plane which is disjoint from [Ao(i?i), 00) and such that fx — 7 
is disjoint from [Ao(i?2), 00), or in other words, such that 7 does not intersect 
[A (-Hi),oo) U (/i - [X (H 2 ), 00)). Suppose also that j(t) = c±t for ±t > T > 
with Imc± > 0. This is illustrated in Figure |l|. The precise values of c± are 
unimportant for our purposes, and indeed there is considerably more leeway than 
this in choosing 7, but the definite separation of 7 from the spectra, ensured by 
Imc± > 0, is crucial. (One can allow slighly subconic separation, but this is of no 
interest here.) Then we claim that 

(2.1) R{n) = -L / fli(A*i)ik(/i - mi) dm. 

To prove this formula, first observe that since 

\\Rj(m)\\ < |Im/x;r\ 



the norm of the integrand in ( |2.l| ) (as a bounded operator on L ) is estimated by 
C(l + |£|) -2 , and h ence the integral converges. Next, note that it suffices to show 



that both sides of (2T) produce the same result when restricted to the range of 
Xi (Hi) <£> Id where I is any compact interval and \i its characteristic function, 
because the union of these ranges is dense. 

Fixing the interval J, the integrand -Ri^i)^^ — M2) is holomorphic for [i\ £ 
IU(ii — spcc(i/2))- Therefore we may deform the contour to one which is the union 
of two curves, 7 and 7, where 7 agrees with 7 for t > T' > and intersects the real 
axis precisely once, somewhere in (sup/, 00), while 7 surrounds / once. Thus 

— / i?i(//i)i? 2 (M - Mi) dfix 

= j Ri{fii)R 2 {n - Mi)^Mi + J i?i(Mi)^2(M - Mi)^Mi- 



(2.2) 



But now, if 7 is moved to infinity, the first integral on the right tends to 0, as 
follows by directly estimating the integral. On the other hand, letting 7 tend to / 
and applying Stone's theorem, the second integral on the right is the same as 



J R 2 {n - Mi) dE x (ni 
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H - spec(H- 2 ) 




H - spec(iJ 2 ) 



Figure 1. Contours of integration in the \i\ plane. 



This is identical to on the range of xi(Hi) 8> Id since applying H — fx to it 

gives (with a slight abuse of notation) 

J^Hx - m + H 2 - (m - Mi))-R2(/x - Mi) dE^fn). 



((Hi - Mi)# 2 (M - Mi) + W) 



jf d£?i(A*i) =X/(ffi)Old, 



i.e. the identity on this subspace. Thus (2.1) is established. 



3. Resolvents of conformally compact manifolds 

In this section we briefly collect some facts about the resolvent family i?M(A), 
which we abbreviate simply as R(X) for the duration of this section, when [M, g) 
is a conformally compact manifold. These are all discussed and proved in |l5[ ], to 
which we refer for all details. 

Recall that M is identified with the interior of a compact smooth manifold with 
boundary M. The compactification M is geometrically natural: M may also be 
identified with the both the geodesic and Martin compactifications, as we discuss 
further below. Locally R(\) is a pseudodifferential operator of order —2, but our 
focus is on understanding the behavior of its Schwartz kernel K(z, z') — K(\; z, z'), 
z, z' G M, when one or both of the variables tend to infinity in M , i.e. to a point of 
dM. This can occur in various ways, of course, and the most efficient way to encode 

2 

this information is to consider if as a distribution on the O-stretched product M 
introduced in jl5| and jy|. This space is obtained from M by blowing up the 
boundary of the diagonal 8Al = d{z = z'}; equivalently, M = [M ;dAi] is the 
disjoint union of M \ 8Al with the interior spherical normal bundle of dAi at 
the corner dM x dM, this set then being endowed with the minimal C°° structure 

2 

containing the lifts of smooth functions on M and polar coordinates around dAi. 

M a has three different boundary hypersurface faces: the left face which covers 
(indeed, is identified with) dM x M, the right face which covers M x dM, and 
the new front face covering dAi. We denote these Bio, -Boi and 5u, and their 
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boundary defining functions pio, poi and pn, respectively. Writing z — (x,y) and 
z' = (x',y') near dM, then 

Pn = \/x 2 + (x 1 ) 2 + \y- y'\ 2 , and p w = x/pu, Poi = x'/pn- 

If dimM = n = k + 1, then write A = £(fc — £), where by convention the region 
Re£ > k/2 corresponds to the resolvent set C \ [fc 2 /4, oo) of A g . Thus for A in the 
resolvent set, 

( = k/2 + i^\-k 2 /4, 
where the branch of the square root is chosen so that its imaginary part is negative. 

Theorem 3.1. (Mazzeo and Melrose, Jl5|, Theorem 7.1]) For £ in the half-plane 
Re£ > k/2, the Schwartz kernel Kq of £)) is a polyhomogeneous distribution 

on Mq with the following properties. First, R(((k — £)) has a decomposition 

R(C(k - 0) = R'(C(k - 0) + R"(C(k - 0), 

where R' is an element in the small calculus (M) of O-pseudodifferential op- 
erators on M and R" is a residual element in the large calculus \I/J~ 00 '''' < '(M) of 
O-pseudodifferential operators. This means that the Schwartz kernel of R' has a 
standard polyhomogeneous singularity corresponding to pseudodifferential order —2 

2 

at the lifted diagonal Alq in M and vanishes to infinite order along Biq and Bqi, 
while the Schwartz kernel of R" takes the form 

p C wP C 01 F", F"er(Mo;4fi ); 

here ir^flo is the lift of the 0-density bundle from the right factor (a non-vanishing 
section of which is given by the Riemannian density dV g ). R{C(k — £)) * s holomor- 
phic, both as a map into the space of bounded operators on L 2 and also into the 
space of distributions on M , for Re C > k/2, and extends meromorphically, as a 

2 

function with values in the space of distributions on M Q , to the complex plane when 
k is even, and to C \ — N when k is odd, with all poles of finite rank. Because it 
is constructed using only the symbol calculus, R' is holomorphic in X. Finally, the 
restrictions of p^K^ to the left face Bio and of p Q i Kq to the right face Bqi are 
nonvanishing for all £ with Re £ > k/2. 

Remark 3.2. It will be convenient below to write 

K C = P C ioPoiF 

where F is smooth on Mq, apart from its conormal singularity at the lifted diagonal 
Ai . 

For us, the main import of this theorem is its conclusion that Kq simple poly- 
homogeneous behavior on the space M . One of our ultimate goals here is to find 
a resolution of the space X 2 , where X is the product of two conformally compact 
manifolds, on which the Schwartz kernel for the resolvent of its Laplacian is also 
simple. 

In the next sections we shall require uniform weighted L 2 estimates for this 
resolvent R(X) as Im A — > oo, and we now show how these follow from the parametrix 
construction in (|l5|]). We set H = — A g here. 
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Theorem 3.3. For any s > and all £ with Re£ > k/2 + s there is a constant C s 
which is independent of Q and such that 

(3.1) R(((k - £)) : x s L 2 (M, dV g ) — ► x s L 2 (M, dV g ) 

is bounded and satisfies 

(3-2) \\RW ~ Q)\W*-v*>») < |ImC ^_ c)| 

for all Q in this half-plane. 

Proof. The proof has a few steps. Setting A = ((k — Q, the parametrix P(X) for 
H — A constructed in pl| satisfies 

P(X)(H — A) = ld+E(X), (H- A)P(A) =Id+F(A), 

where -E(A) and -F(A) are residual, and is related to -R(A) by the formula 

(3.3) R(X) = P(X) - E(X)P(X) + E(X)R(X)F(X). 

The desired uniform estimate for R(X) then follows from the standard uniform L 2 
estimate 

II^A)!!^,^) < — — , AeC\spec(i7), 

and from appropriate uniform estimates for -P(A), E(X) and F(X) in weighted spaces 
which we now derive. 

The main work involves demonstrating a uniform estimate in weighted spaces 
for the resolvent of the Laplacian in hyperbolic space, and so we turn to this first. 
Uniform boundedness of 

R(C(k - C)) : x s L 2 (U k+1 ; dV) — > x s L 2 (U k+1 ; dV) 

is equivalent to the uniform boundedness of the conjugated operator x~ s R(£(k — 
C))x s on L 2 (U k+1 ;dV). The Schwartz kernel of this conjugated operator is 

K s c (z,z') =K a c (z,z')(x'/x) s =K c °(z,z')(PoiMo) s , 

where K®(z,z') is the Schwartz kernel of the resolvent of the Laplacian in H fe+1 . 
Recall from ]l5| that there is an explicit formula for K®: 

(3.4) K$(5) = 

k-2 

1 d \ 2 ( 1 _-(£-fc/2),5 



sinhi5 9(5/ \sinhi5 /' 
e -(C-k/2)w( CQshuj _ cosh( 5)- fe / 2 dujj k odd. 



Cfc 

Here we use that the resolvent is a point-pair invariant, so R(X; z, z') only involves 
the Riemannian (hyperbolic) distance between z and z 1 , and equivalently, may be 
written in terms of the elementary point-pair invariant S(z, z') defined by 

\z - z'\ 2 

cosh5(z,z') = 1+1 -L. 

2xx' 

There are two regions in (H fc+1 )Q where behaves slightly differently as £ — > oo. 
These are one near the diagonal, say 8(z, z') < 1, and one away from the diagonal, 
say (5(z, z') > 1. We use a partition of unity to divide the Schwartz kernel into 
two parts. It suffices to prove L 2 boundedness of each piece separately. In the 
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former region, S(z,z') < 1, the factor (poi/y°io) s is boun ded, independent of £, 
hence the uniform L 2 boundedness of K^, and thus (3.2) can be proved just as 



the L 2 boundedness of K®, directly from properties of the Schwartz kernel. More 
explicitly, the result in this region may be deduced either by invoking the standard 
extension of the symbol calculus with spectral parameter, or else simply by direct 
calculation. On the other hand, in 5(z, z') > 1, P| decays exponentially as £ — > oo, 

so the additional factor (poi/pw) s , which takes the form e s ^^ p ' z ^ s ^ p ' z where p 
is some fixed point in the interior of H fc+1 , does not make much difference: the 



Schwartz lemma shows the desired bound (3.2) (and in fact better bounds for this 
piece!). 

To proceed, we recall that P(A) is constructed in stages, and as a sum of two 
terms, P(A) = Pi (A) + P 2 (A). Here Pi (A) is in the small calculus, and is a slight 
modification of the operator P[(X) obtained from the use of the symbol calculus 
to solve away the conormal singularity along Alq. By standard methods, such 
a P{(A) may be constructed so as to depend holomorphically on ( G C, and to 
have uniformly bounded norm on L 2 . It also acts on x a L 2 for each s, with norm 
depending on s, but not Q. Next, p2(A) is obtained after solving away the first I 
terms of the Taylor series expansion of the error term Pi (A) = I — (H — A)P 1 '(A), 
Pi (A) G °°(M), where I is sufficiently large (greater than s). This approximate 
solution of (P — A)P 2 (A) = Pi (A) is found as follows: first solve the normal equation 

(N(H) - X)N(P 2 (X)) = N(F 1 (X)), 

which is the restriction of the equation to the front face B\\. But N(H) is naturally 
identified with the Laplacian on H fc+1 , and so we may apply the result of the 
discussion above to choose P 2 (A) with this property, which satisfies the estimate 
and such that P 2 (A) = Pi (A) — (H - A)P{(A) vanishes to first order along 
B n . Thus 

(H \){P[{\) + P^(A)) = / - P^(A), P^(A) g *-°°'^(M). 
Applying the first t terms of the Neumann series for (/ — P^A)) -1 , 

/ + S w = 1 + P 2 '(A) + . . . + (P^A))^ 1 
on the right on both sides of this equation yields 

(PJ-A)(Pi(A)+P 2 (A)) = /-P 2 (A), 

where 

P j (X)=P;(X)(I + S^(\)), 3 = 1,2, and P 2 (A) G *„ °°' C+ ^(M). 



Clearly 



np 2 (c(fc-c))ii BW 3)< | ImC( fc_or 



with C independent of £. We also obtain a remainder term P(A) from applying 
H — A on the right of P( A) = Pi (A) +P 2 (A) , and this satisfies uniform bounds (from 
L 2 to x s L 2 ) of the same form. 

These facts taken all together finish the proof. □ 
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4. Analytic continuation of the resolvent of H 



In this section we shall show how to use the integral formula (2.1) to obtain an 
analytic continuation for the resolvent Rx(fi) = {H — fi) on X = Mi x M% past 
the continuous spectrum of H = — A g = H\ + H 2 - For this continuation we can 
either regard Rx (n) as an analytic family of bounded operators between weighted 
L 2 spaces, or else view its Schwartz kernel as an analytic function with values in 
an appropriate space of distributions. The key ingredient here is the existence of 
similar analytic continuations for the resolvents Rj{fi) — {Hj — /it) . Although this 
continuation result holds in considerably greater generality than just for products 
of conformally compact spaces, we shall focus exclusively on this case for the sake 
of being specific. 

We first set up some notation. Let dim Mj — kj + 1. Then it follows from 
Section cf. also Q that the spectrum of Hj decomposes into the union of a 
band of continuous spectrum [kj/4, +oo), as well as possibly a finite number of L 2 
eigenvalues Xji, i = 1, . . . , Nj, j = 1, 2, in (0, k 2 /A), with the corresponding finite 
rank eigenprojections Hji- Next, each Rj(fij) continues meromorphically from the 

resolvent set to the Riemann surface E for yj fij — k 2 /4, which we think of as two 

copies of C attached in the usual way along a cut Cj extending from k 2 /A to oo; the 

resolvent set is identified with the subset of E where Im J jij — fc|/4 < 0. Usually 

Cj is taken to be the ray along the positive real axis, but it will be convenient to 
choose it differently later. As already noted, this continuation of Rj is either as a 
map into an appropriate spaces of distributions, or else for any given s > 0, in the 

region Im [ij — k 2 /4 < s, as bounded operators from x s L 2 to x~ s L 2 . In any case, 

its poles are of finite rank, and we denote them by Xji, with corresponding finite 
rank residues EL^ (the tildes are meant to distinguish these from the eigendata for 
Hj). Note, however, that these poles of Rj in the nonphysical part of E need not 
be simple. 

We assume until near the end of the argument that neither Hj has any L 2 
eigenvalues. Define k by 

k 2 h 2 h 2 

T = T + T' 

Then R(fi) is analytic in C \ [k 2 /4, oo), and we wish to show that it continues 
analytically past the continuous spectrum. We do this by deforming the contour of 
integration 7 in (|2.l| ) in the following manner. First fix in the resolvent set, so that 
Ri(pi)R2(p — fJ-i) is defined and analytic for fix ^ C\ U (fj, — C2), i.e. outside of two 
horizontal rays, one extending from kf/A to the right and the other from fj, — fcf/4 
to the left. Next, rotate these cuts to rays Cj by pivoting them by some angle a 
counterclockwise around their endpoints. Thus we have 'exposed' two sectors of 
angle a from the nonphysical portion of E, and at the same time concealed an 
equal portion of the physical part of E. Now it is possible to deform 7 to a new 
contour 7' which lies partly in the newly uncovered sectors in C \ (C\ U C2), as in 
Figure |^. Finally, /1 can then be moved into the nonphysical region. 

We make this a bit more explicit. Suppose that /x is such that /x° — k\/A ^ Ah 
and (U° — k 2 /A ^ A2;' for any i, i' , and also fi° ^ Xu + X 2 i> ■ Rotate the cuts Cj by 
an angle a > arg(^° — k 2 /A), j = 1, 2. Now deform the contour 7 past fi° — k\jA to 
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Figure 2 . Contour of integration to describe the analytic contin- 
uation of R(n). In the second picture, the cuts C\ and C 2l cor- 
responding to the spectra of H\ and i? 2 , respectively, have been 
rotated to the cuts C\ and C 2 . Thus near the point /u? labeled here, 
Ri(Hi) is not a bounded operator on L 2 (Mi). Now fi, and hence 
C 2 as well, can be moved upwards in the direction of the arrow. 



a new contour 7', such that this point is now below 7'. During the deformation, 7 
may pass through a finite number of poles of R\, yielding residues nii®.R 2 (/x— Ah). 
So long as fi is (sufficiently) far away from C 2 , fi — fJb\ is in the resolvent set of 7? 2 , 
hence no poles of R2 are encountered during this deformation. Thus we have a new 
representation 

r (h) = 77— I Ri{ni)R2{p- - mi) d m - V" flu <g> i? 2 (M - Ah). 

2771 J* fit^e 

Since fi is in the resolvent set, the left hand side is still a bounded operator on L 2 , 
although on the right hand side, one factor of the integrand Ri(fii) is not bounded 
on L 2 along the entire contour. 

At this point we merely have a new representation of an operator we already 
know exists, namely R(/j,) for fi in the resolvent set, in terms of operators which 
are only bounded between weighted spaces. However, fixing 7', we can now let [i 
vary arbitrarily below this curve, and hence across the spectrum [fc 2 /4, 00), as long 
as the ramification point of R 2 does not lie on 7'. When fj, — A 2 i crosses 7', 
the residue term — i?i(/x — A 2 i) ® Uu is produced. This yields 

R{fj) = J Ri{n\)R2{n - Hi) dfi\ 
- ^2 fl u <g> R 2 (n - Ah) - J! ~ ^) n 2i ; 

finite finite 

Each of the residue terms here continue meromorphically in /i past the spectrum, 
with ramification points at Ah + fc 2 /4 (since R2 is ramified at fcf/4) and similarly, 
with the indices 1 and 2 interchanged. These terms also have poles at Ah + A 2 j, 
with finite rank residues. 
We have proved 
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Theorem 4.1. The resolvent R(p) for H on X = M\ x M2 extends across the 
cut [k 2 /4,oo) to a meromorphic function (with values in an appropriate space of 
distributions) on a Riemann surface T ramified at Xu + fcf/4 an d ^2% + (these 
points are known as Regge poles), and with finite rank poles at Xu + X%i- 

When either Hi or H2 has L 2 eigenvalues, then the only difference is that 7' 
must cross these as well. Thus, the analytic continuation can be written as above, 
but now we must include sums over the Xji too, and there are new ramification 
points at Ah + fc|/4 and X 2i + fcf /4. 

5. COMPACTIFICATION CONSTRUCTIONS 

We now turn to the other major theme of this paper, which is the various ways 
one might compactify the conformally compact spaces Mj or their product X — 
All x M2. The general problem of finding good compactifications of Riemannian 
manifolds, and in particular of (locally) symmetric spaces, has been an area of 
active research. We refer to JTT| and Q| for more discussion of this in the symmetric 
space setting, and to and [|22| and [[| for some beautiful results in some general 
geometric settings. 

There are two different compactication constructions we shall discuss here, the 
geodesic compactification (sometimes also called the conic compactification) , as 
well as the Martin compactification. Each has a key role in understanding different 
aspects of the global geometry and function theory of a space. We define these now 
in turn. 

When Z is a complete, simply connected manifold of nonpositive curvature (i.e. a 
Cartan-Hadamard manifold), then the geodesic compactification Z is obtained by 
adjoining to Z an ideal boundary dZ, points of which arc equivalence classes of 
geodesic rays. Two geodesic rays ^(t), j{t), t > 0, are said to be equivalent if 
d(j(t),j(t)) remains bounded as t — > 00. Thus in R™ any two parallel lines are 
identified, as are any two geodesies in the ball model of hyperbolic space H™ which 
converge to the same point on If q G dZ, then a neighborhood system at q is 

given by sets of the form {Z\Br{p')) n exp p ([T, 00) x U), where U is an open set in 
the unit sphere S p Z in T p Z. A point on this set is on some geodesic j(t) emanating 
from p with 7'(0) € U and t > T as well as in the exterior of the ball Br(p'). 
Thus whenever Z is Cartan-Hadamard, then Z is homeomorphic to a closed ball 
B n , and its boundary dZ is identified with the unit tangent sphere S p Z for any 
p G Z . There are fairly obvious modifications of this construction when Z is not 
necessarily simply connected, or only has nonpositive curvature outside a compact 
set, and then the structure of Z is more complicated, but only on account of the 
topology of int Z. 

It is of substantial interest to understand when the compactification Z carries 
more structure than its initial definition as a topological space. In particular, we 
would like to determine when Z is naturally defined as a smooth manifold with 
boundary, or with corners. For example, When Z = R" or H™, then as we have 
already stated Z w B n , but obviously in these two examples, Z is 'really' a smooth 
closed ball. An interesting and underappreciated feature of this construction is that 
although it is possible to identify each of these compactifications with a ball, the 
smooth structures at the boundary are different, as we now explain. 

When considering whether Z admits a natural smooth structure, the first key 
point is the regularity of the transition maps, as we now describe. Remaining within 
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the context of Cartan-Hadamard manifolds for simplicity, the transition maps are 
the homeomorphisms between the unit tangent spheres at any two points p,p' G Z, 

SpZ > dZ < Sp'Z 

provided through the geodesic spray from these two points. When Z has curvatures 
bounded between two negative constants —a 2 and — b 2 , then from these tran- 
sition maps are of Holder class C 0,a , a = a/b, and accordingly in this generality, 
Z has only a Holder structure. However, for either of the examples above these 
transition maps are smooth; this is also true for general conformally compact man- 
ifolds for suitable localized versions of these transition maps (see The other 
key point concerns the choice of a class of defining functions for dZ (or, if Z is to 
be regarded as a smooth manifold with corners, then for subsets of it which are 
identified with the various boundary hypersurfaces) . Smoothness of the transition 
maps and choice of defining functions together determine the C°° structure on Z . 

There is considerable flexibility in choosing (equivalence classes of) defining func- 
tions. For W 1 it is most natural to use the radial compactification with defining 
function for <9R™, pi = l/\z\ = l/dist(z,0), z € M™, the inverse of the polar dis- 
tance variable. On the other hand, the Poincare ball model for H™ suggests that the 
more natural choice now is p 2 = exp(— dist (z,0)), any fixed point in HP. These 
two defining functions for dB n are quite different, since p\ = — l/logp2- We say 
that (B n ,pi) is the 'log blow-up' of (B n ,p 2 ). For later reference, we shall refer to 
defining functions defined as the reciprocals of the Riemannian distance function 
or its exponential as being of polynomial or exponential type, respectively. 

The difference between these polynomial and exponential type defining functions 
appears most clearly in the spaces of functions with polyhomogeneous behaviour at 
the boundary, since this involves expansions in a discrete set of complex powers of 
the defining function p along with nonnegative integer powers of its logarithm, log p. 
In particular, functions which are polyhomogeneous with respect to exponential- 
type defining functions, with all exponents having positive real part, are residual, i.e. 
rapidly decreasing, with respect to the log blow-up structure. The guiding principle 
for which defining functions to choose as primary is determined by the asymptotic 
behaviour of eigenfunctions of the Laplacian; one prefers eigenfunctions not to be 
automatically residual! The choices made above for 1™ and H™, respectively, are 
vindicated by the analysis of manifolds with 'scattering metrics', see and with 
conformally compact metrics, see |1| and | fl3|| . 

As indicated above, the geodesic compactification of the conformally compact 
manifold (M, g) is just M with its usual smooth structure. It is obviously of interest 
to identify the geodesic compactification for the product of two such manifolds, 
X = Mi x M 2 . In fact, dX is the simplicial join of dM 1 and dM 2 . More specifically, 
it is obtained from 3M\ x dM 2 x [0, n/2]g by collapsing each dM\ x {q 2 } x {0} 
and {qi} x dM 2 x {ir/2} to a point. Alternatively, it can be described as the 
manifold with corners obtained by blowing up the corner dM\ x dM 2 in the product 
Mi x M 2 , and then blowing down the M 2 fibres in 8M\ x M 2 and the Mi fibres in 
Mi x 8M 2 in the original faces. There are three boundary hypersurfaces, dM\ x M 2 
and Mi x dM 2 , corresponding to limits of geodesies of the form ("fi(t),p 2 ) and 
{pii J2(t)), respectively, and also the new face dM 1 x dM 2 x (0, ir/2) corresponding 
to limits of geodesies (71 (at), 72 ((3t)), a 2 + 1 — 1. It is less clear what defining 
functions to use for these faces, though since the 'new face' covering the corner 
comes from the two-dimensional flats, i.e. products of geodesies in each factor, it 



14 



RAFE MAZZEO AND ANDRAS VASY 



is reasonable that we should use a polynomial type defining function here. In fact, 
using the usual coordinates on each factor, we shall use 



In other words, the smooth structure on X is the one induced from the normal 
blow-up at the corner of the log blow-ups of the two factors Mi and Mi- The 
reason we use the log blow-up at all faces is that as we shall prove in the next 
section, the asymptotic behaviour of the resolvent involves expansions in powers of 
these particular defining functions. 

The other compactification we consider here is due to Martin , and uses the 
function theory of A g to associate a set of ideal boundary points dZ to Z . Notably, 
it may be carried out in great generality for pairs (Z, H) where H is a semibounded 
self-adjoint elliptic operator on a space Z, though we shall always assume here that 
H is the Laplacian. Let Ao = inf spec(iJ). Then for every A £ M \ (A , oo), there is 
a compactification Zm(A). Actually, it follows from the construction that 2m(A) 
is identified with Zm(A') for any two numbers A, A' < Ao, so one needs to consider 
only Zm(^o) and any other Z^(A). 

Fix A < Ao- By [g3|, the set of solutions u to the equation {H — X)u = which 
remain everywhere positive is nonempty; we denote this K + -invariant set by V+(X). 
The structure of this positive cone is encoded in the slice 7- + (A) = {u £ 7-+ (A) : 
u(p) = 1 for any fixed p G Z. It follows readily from the Harnack inequality and 
elliptic estimates that for any sequence Uj £ 7- + (A) there is a subsequence Uj> 
converging to a point in this space, and so 7- + (A) is compact. It is also obviously 
convex. Let £ denote the set of its extreme points. Then for any u S V^.(X), 
the Krein-Milman theorem gives a measure dm u (e) supported on £ such that u = 
J £ dm u (e). This is the generalized Poisson representation theorem! The set £, or 
£ (Z, A) is called the minimal Martin boundary of Z. If u £ £, then whenever 
v € Va_ (A) and v < u then v — cu for some constant < c < 1, and this justifies 
the moniker 'minimal'. 

If Z were to naturally embed in 7-^ (A), then the closure of its image would be 
an obvious way to compactify it. Unfortunately, this is not the case, but instead we 
consider the resolvent kernel i?(A; z, w). This is a solution of (H z — A)i?(A; z, w) = 
when z ^ w, but is singular at z = w and in addition, R(X;p,w) =/= 1. Thus we 
define 



so that u w (jp) = 1 and u w is a regular solution of {H — X)u w = on X \ {w}. Now 
let Wj be any sequence of points in Z which leaves any compact set. Then some 
subsequence u w >. converges to an element of 7- + (A). The (full) Martin boundary 
is the set of equivalence classes of these sequences, or equivalently, is the set of 
all possible functions u obtained as limits in this fashion. We label the different 
boundary points by q G d]tfZ(X), and write u q {z) for the corresponding limiting 
solution. The Martin compactification Zm (or Zm(A)) is the union of Z and 
8mZ{\). There is a metric on the set of functions u w {z), w G Zm, given by 



Pj 



1/logXj, j = 1,2, 



/■ = 



\{P1,P2)\- 



u w (z) — i?(A, z,w)/R(X,p,w) 




Thus Z m not only a topological space, but a metric space. 
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This definition must be modified when A = Ao and there is an L 2 eigenfunction 
u with eigenvalue Ao, since then the resolvent kernel R(\q,z, w) does not exist. 
In this case u > 0, and in fact R+ • u = P+(Ao), i.e. the only positive solutions 
are positive multiples of uq. Hence it is consistent to let Zm(Ao) be the one point 
compactification of Z then. Otherwise, if Ao is not in the point spectrum, then the 
definition is the same as before. 

If Z = R™, then Ao = and it is well-known that c?m R™(0) is a single point, so 
that (R™)m(0) is the one-point compactification S n . On the other hand, for A < 0, 
the extreme points of V+(X) are the exponentials e x |£| 2 = —A, and this sphere 
of radius V^A is the full Martin boundary, and so (1") M (A) = B" . If Z = H" 
and n = k + 1, then Ao = fc 2 /4 and Zm(A) is the closed ball B n for all A < Ao. 
The minimal positive eigenfunctions are the ones of the form x k / 2 +\f k ~ / A ~ x f or a n 
possible choices of upper half-space coordinates (i.e. choice of which point on the 
boundary of the ball to send to infinity) . Fro m Jjl| it follows that that when M is 
conformally compact, Mm (A) is still equal to M, To survey other cases relevant to 
us in which the Martin compactification is known, when Z is Cartan-Hadamard with 
curvatures pinched between two negative constants —a 2 and — b 2 , then Anderson 
and Schoen |2| and Ancona [Q proved that Z M (A) = B n (as metric spaces) . There 
has been recent significant progress in determining the Martin compactifications 
for general symmetric spaces of noncompact type; definitive results are proved in 
the recent monograph |7j, and this has an extensive bibliography of the literature 
on these developments. 

Of particular relevance to us here is the work of Giulini and Woess || where 
the Martin compactification of the product of two hyperbolic spaces is determined. 
Their result is that now A M (A) is identified with p^xlf^; dW^ x 9H™ 2 ] , the blow- 
up of the product of the two balls along the corner; recall that this space appeared 
as an intermediate picture in the description of the geodesic compactification. Their 
proof uses rather involved global heat-kernel estimates, and one of the motivations 
of this paper was to demonstrate how this result, and the analogous one for products 
of conformally compact spaces, may be obtained in a more straightforward manner 
using resolvent estimates and stationary phase. 

To conclude this section on compactifications, we state once again that our pri- 
mary interest is in obtaining a compactification of the double-space A x A, where 
X = All x M2, which is natural with respect to the resolvent. More specifically, we 
wish that at least for A in the resolvent set, -Rx(A) should have at most polyhomo- 
geneous singularities at the boundary hypersurfaces of this compactification (apart 
from its usual diagonal singularity) . This compactification is determined by exam- 
ining the asymptotics of i?v(A; z,w) as the points 2 = (2:1,2:2) and w = [w\,W2) 
diverge in all possible directions. We begin the study these asymptotics in the next 
section. The Martin compactification is, at best, essentially a 'slice' of this resol- 
vent compactification, and in any case is easily determined from this asymptotic 
analysis, as we do in the final section. We shall see there that when either Hi or H2 
has eigenvalues below the continuous spectrum, then the Martin compactification 
of A is obtained by substantially collapsing part of the boundary of the geodesic 
compactification A, and so loses a lot of information about the fine structure of the 
resolvent. Thus the resolvent compactification is the primary object of interest. 
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6. ASYMPTOTICS 

In the remainder of this paper we shall give a more detailed description of the 
structure of the resolvent on X. As a first approach to this we adopt the 

more traditional viewpoint and derive the asymptotic behavior of when 
/ G C°°(X). (Functions vanishing to all orders at dX are a suitable analogue of the 
space of Schwartz functions.) More generally, it makes absolutely no difference if we 
allow / to be the sum of an element of C°°(X) and a distribution of compact support. 
In particular, all of the calculations below apply when / = S p , p G X = Mi x Mi\ 
indeed, this is the basis for our identification of the Martin boundary of X. However, 
we shall simply assume that / is Schwartz, and also, in this section, that /i is in 
the resolvent set for H. 

Recall our convention that when \ij is in the resolvent set for Hj, then the 



imaginary part of J fij — kj /4 is negative, and by the results of last section, Rj (fij)f 
is then of the form 

x j v .</;• 9 J G C (Mj). 

Also, if cfiji is an I? eigenfunction of of Hj with eigenvalue Xji, then 



<i>ji =x- v gj, fi, € C [U j): 

the Schwartz kernel of Hji is <f>ji(zj) <g> 4>ji(z'j) if the eigenvalue is simple, and is a 
finite sum of such terms otherwise. 

We start with a slightly stronger result concerning the structure of i£i2(Mi> M2) = 
Ri(^i) <S> i?2(M2) applied to / e C°°(X), considered as a function of both fii and 
/i2 . Although this follows directly from the corresponding statement in each factor 
if / = /l ® /2j fj G C°°(Mj), it is convenient to prove the general result directly. 
Thus, the kernel of i?i2(A < i, M2) is polyhomogeneous on (Mi)g x {M-^'q- Let ttl, 
7Tfl be the projections to the left and right factors of Mi x M%. That is, if ttlj, 
resp. 7Trj, denote the projection of (Mj)q to its left factor, resp. right, factor, then 
til = ttl.i x 7Ti 2, tt_r = ^ R,i x ^r,2- Note that itl, ttr, are b-fibrations. Then 
R\2{n\, M2) applied to / is given by the push-forward of Ri2{^\, ^2)71"^/ under the 
map ttl- Thus 



(6.1) #12(^1,^2)/ = a^i a: 2 v 9,9^C (Mi X M2). 

The dependence on and /i 2 her e is uniform in the strong sense that the function 



g(zi, Z2, A*i, M2) appearing in (6T) satisfies 

.9 e C 00 ^ x M 2 x (C \ spcc(i? 1 )) All x (C \ spec(ff 2 )) M2 ), 

with natural extensions corresponding to the meromorphic extension of the 
This follows directly from the usual push-forward formula Although we use 

this result throughout the section, we usually state arguments for simplicity as if 
Ri and i? 2 are applied separately to f . In one particular case, when analyzing 
R{ti)f for /1 real, below spec(-ff), and when either Hi or H2 have L 2 eigenvalues, 
we need a stronger result, where / is not required to be Schwartz. We postpone 



that discussion until it is required, see the arguments prececding Theorem 6.3. 

The asymptotic behavior of R(fj,)f(z) must be analyzed in three separate regions: 
near dMi x M2, near Mi x dM2 and near the corner dM 1 x dM2- Using coordinates 
z — (zi, Z2), Zj — (xj, y.j), these correspond to x\ — > 0, X2 > c > 0, or xi > c > 0, 
X2 —> or xi,X2 — * 0, respectively. 
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6.1. Asymptotics at Mi x dMi and at dM\ x M 2 . We first describe the uniform 
behaviour of R(n)f on Mi x M2, i.e. at infinity in M2. This is given as an asymp- 
totic series in powers of — 1/ log x 2 ; we shall think of this later as an asymptotic 
expansion on the logarithmic blow-up of M2, which simply means that we change 
the C°° structure of M 2 by replacing the defining function x 2 for the boundary by 
-l/logx 2 . 

Suppose first that H\ has no L 2 eigenvalues. We analyze using the 

representation fl2.l| ) for R(p), by shifting 7 so that it passes through inf spec(i?i) = 
fci/4, and so that the minimum of Im y/fi — k 2 /A — \x\ along this path is attained 
at that point, see Figure ||, and then applying (complex) stationary phase. 




Figure 3. Contour of integration, 7, used to describe behavior in 
Mi x M2 . Here c = Im \j [i — k 2 /4, and u 2 = yjji— — Ml- 

To set things up for stationary phase, first note that R 2 {^, 2 ) maps C°°(M2) to 
a .^/ 2 +Vw- fc i/ 4 C oo(]g 2 )_ Therefore, the oscillatory factor e (W2+V^- fe i/ 4 ) lo s^, 
l-i 2 = [i — fJ-i appears in the integrand. On the other hand, Ri(pi) is not smooth 
at /ii = &i/4, but instead has an expansion in powers of \J fii — fcj/4. Thus it 
decomposes into the sum of odd and even powers, respectively, Ri(fii) — R° dd ([ii)+ 
R1 vcn (m), and we write Rf d (^i) = y/JTi - kJ/AG^fii-kf/A) where d is smooth. 
Because the phase y/fx — (Ti — k 2 /A logx 2 is not stationary at fci/4, the smooth part 
i?° vcn contributes only terms decaying faster than any power of — 1/ log %2- Thus 
we are left only with 



(/ E2 /2 + VM-w-fel/4)io ga y /ii _ fc f/4Gl(/il - fc?/4)^i,Z,/)d/ii, 



e 

7 

where g is smooth in all variables. Actually, it suffices to take this integral only over 
some compact segment of 7 containing fef/4 and where the decomposition of R\ 
into even and odd parts is valid; using the uniform weighted L 2 estimates of §3, the 
integral over the remaining portion of 7 contributes a term vanishing like x 2 2 ^ 2+m , 
m > can be made as large as wished by increasing the length of the compact 
segment, and is therefore negligible in the asymptotics. Changing the variable of 
integration to r = yjx\ — k 2 /4, so that /ii = k\/A- 



e (/= 2 / 2+ V^W^)logx 22r 2 G(r2) ~ (fc 2 /4 + T 2 ; Zj zl) dT _ 
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The phase function 4>{t) = wJjl — fc 2 /4 — r 2 log x 2 is now stationary at r = 0; 
furthermore, </>"(0) = —1/Wfi — k 2 /4, and so the further change of variables a = 
t^J— log x 2 j (/i — fc 2 /4) 1 ' 4 reduces this integral to 



I do- 
where g is obtained from g by replacing a by t. Applying the stationary phase 
for phase functions with nonnegative imaginary part, see p0[ Theorem 7.7.5], we 
obtain that 



(6.2) R(»)f = a ; 2 fe2/2+^ ^ /7 ^(-l/logx2) 3 / 2 ff , 

where g is C°° on the logarithmic blow-up Mi x (M2)i og of Mi x M2, i.e. in the 
variables (21, —1/ log X2, 2/2) for xi > c > 0. (Of course, g is continuous on Mi x M2; 
it is the lower order terms in the asymptotics that necessitate the change of the 
smooth structure.) In fact, 



9\m M = C (/,-fc 2 /4)3/V 2/2 ^^^G 1 (0)i? 2 (M-fc?/4)/ 



MiXdM 2 



where c is a constant arising from the stationary phase lemma. 

There is a final simplification of this formula, arising from the fact that we can 
identify the operators here slightly more explicitly. First, according to G3], for each 
zi € Mi, 



: k2/2 '^^ k2/i R 2 ( P -k 2 2 /A)f 

P*( M - fcf/4; y 2 , 4)/(*i, 4) ^2 = ^(/* " ^/ -I )(/!(: i . </j ) • 

fcf/4) is the Poisson transform on A 
H — fcf/4, and Pf its transpose. Next, by definition, 



ldM 2 

where P 2 ([i — is the Poisson transform on M 2 at the spectral parameter 



i?!(fcf/4 + T 2 ). 



We can already see from how it was introduced that this operator arises from 
the non-holomorphic part of R\, and so its Schwartz kernel is smoo th on (Mi)q 



because the diagonal singularity of R\ , represented by R[ in Theorem 3A , does not 



contribute to it. However, we can see this more directly by differentiating 

(H 1 -k 2 1 /4-T 2 )R 1 (k 2 1 /4 + T 2 ) = I 

with respect to r and setting r = to get (Hi — fc 2 /4)G(0) = 0, and so concluding 
that Gi(0) is smooth by elliptic regularity. The Schwartz kernel Gi(0; Z\, z[) is a 
familiar function on hyperbolic space: it is a non-square-integrable eigenfunction 
for H\, with threshold eigenvalue k 2 /A, invariant under the group of rotations which 
fix z'i, and is known as a spherical function. We shall denote it by <Si(fc 2 /4; z\, z[), 
and regard it as a Schwartz kernel. Altogether, we have now shown that 

(6-3) 5lM lX9 M 2 =c(M-fc 2 /4) 3/4 5i(fcf/4)Pi(M-fcI/4)/- 

When Hi has L 2 eigenvalues, then shifting the contour 7 through these gives a 
contribution — Ilii ®R 2 ([i — An) from the residues of the integrand. (For simplicity 
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we assume that the ramification point k\ /4 is not a pole of the meromorphic con- 
tinuation of R\ to the Riemann surface, though this can obviously be handled too 
by the arguments preceeding Theorem [T]^.) Hence in this case we obtain 

(6.4) 

D , w k 2 /2+i^f n-k 2 /A , , .3/2 

«(m)/ = x 2 v {-l/\ogx 2 y' g 

Nl ' 

+ ^4 a/2+i V"- A "-^ /4 fli) where g, 9i € C°°(Mi x (M 2 ) log ). 

In fact, gi = (f>u ® g^, with <?j G C°°(Af2) and (f>u an eigenfunction of ifi with 
eigenvalue A H ; indeed, ^| MlX gM 2 = ® P 2 *(m - Ah))/- 

In this region, near Mi x <9M2, only the L 2 eigenvalues of Hi play a role in the 
asymptotics, but perhaps surprisingly, not those of H 2 . Notice that all terms in 
( |6.4| ) which come from the L 2 eigenvalues of Hi dominate the term coming from 
the continuous spectrum since Im \JJi — W/i < Im \J [i — \u — k\jA. In addition, 
the term corresponding to the lowest eigenvalue A10 of Hi dominates all the other 
terms. This will be important later in the determination of the Martin boundary 
in the presence of bound states. 

The asymptotics of R(n)f at the other face dMi x M 2 is completely analogous. 
The same calculations lead to the expression 

(6.5) 

R[n)f = x x v (-l/logxi) J/ g 

N\ . 

+ g^ l/2+ VM-A 2i - fc ?/4^ where g g coo^^ x Ma) _ 

i=l 

In addition, 

(6-6) 9\dM lxM2 = c(M - k 2 /4) 3 ^P^ - kl/A)S 2 {kl/A)f, 

where P[ is the transpose of the Poisson transform on Mi and 5*2 is the 'spherical 
function' at eigenvalue k 2 /A for M 2 . 

We note once again that even if both Hi and H 2 have L 2 eigenvalues, only those 
of H 2 contribute to the asymptotics when xi — > 0, x 2 > c > and only those of 
Hi contribute to the asymptotics when x 2 — * 0, xi > c > 0. This indicates that 
the asymptotics at the corner dM\ x dM 2 , where both x%,X2 — ► must be more 
complicated, at least in the presence of bound states, because it intermediates this 
transition. 

6.2. Asymptotics at dM\ x dM 2 in the absence of L 2 eigenvalues. We now 

proceed to the analysis of R(fi)f at this corner. We have already seen that necessity 
of logarithmic blow ups at the side faces of the product. Correspondingly, the 
asymptotics at the corner necessitate that we pass to the log blow-ups of both 
factors Mj, which we denote (Mj)\ og . In terms of these, define 

X = [(Mi)iog x (M 2 )io g ;9(Mi) log x 9(M 2 )io g ]. 

This is the correct space on which to consider asymptotics of R(fj,)f, as we shall 
now prove. 
We let 

Pj = -1/logatj-, j = 1,2, 
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be boundary denning functions of (Mj)i og , and we keep denoting their pull-back to 
X with the same notation. Also let yi, 1/2 denote coordinates on d(Mj)i og , j = 1, 2. 
Thus, X is the blow-up of (Mi)i og x (M2)i og at the corner pi = P2 = 0. Hence 
valid coordinates in the region where P1/P2 < C, for some C > 0, are given by 

P2, s = pi/p 2 , Vi, and y 2 . 

In terms of the original boundary defining functions Xj this means that in any 
region where loga^/ logxi < C, for some C > 0, we use the projective coordinate 

-l/logxi logx 2 

s = = 

-l/logx 2 logxi 

along the new 'front' face of X covering the corner, and — l/logX2 as a defining 
function for this face. Thus s — > upon approach to the lift of dMi x M2, while 
s — > 00 on approach to the lift of Mi x DM2- We note that a total boundary 
defining function of X is given by 

I -2 1 -2\— 1/2 
P=KPl +P2 ) . 



i.e. with the usual Euclidean notation, rj = Pj , r = p 1 , r = + r| . This 
explains the appearance of p -1 in our results below. In the region s < C, p\ = 
— 1/ logxi is another total boundary defining function; it will be used for the initial 
local calculation and then we change to p for the global statements. The trivial 
identity 

x 2 =e losX2 = e slosXl =x{ 

will be used repeatedly. 

The integrand Ri{p.\)R2{p-~ Pi) f in the contour integral representation of R(p)f 
has the form 

x 1 x 2 g, 



or equivalently, 



(6.7) exp 



ki/2 + sk 2 /2- 



logxi \g 



where g is C°° in p, pi and on Mi x M%. The expressions inside these square roots 
assume values in C\ [0, 00), and we assume the square roots have negative imaginary 
parts in this region. 

Suppose first that neither Hi nor H2 has L 2 eigenvalues. We again do a station- 



ary phase analysis of the integral of (6.7), and for this it is necessary to choose the 



contour of integration so that, if we set 

F(p 1 )=F s {p 1 ) - yjm- kl/A + s^p-pi~ kl/A, 

then the supremum of ImF(pi) along 7 is as negative as possible. Both p and s 
are parameters here, and we must choose the contour differently corresponding to 
the different points of the front face of X, i.e. the different values of s. Since F is 
analytic outside [fcf/4, +00) U (p — [k^/A, +00)), its critical points and those of its 
imaginary part are the same. In fact, F has a unique critical point, located at 

(6.8) rf-rfW- "-^,''^ 4 . 

1 + s z 
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and this is a saddle point of the harmonic function Ini-F^i). (We return to this 
and shall explain it in greater detail below.) Moreover, this critical point always 
lies on the straight line segment connecting the two ramification points fc 2 /4 and 
[i — fc|/4, and tends to the former as s — > +00 and to the latter as s — > 0. Next, 

F04) = V(M-fc74)(i + s 2 ), F'( M ?) = 0, 

and 

F"( M °) = - i S - 2 (l + s 2 f/ 2 {ii - fc 2 /4)- 3/2 . 

Therefore we may choose an appropriate contour 7 so that IvnF(^(t)) attains a 
non-degenerate maximum when j(t) = The stationary phase lemma may now 
be applied as before, although the singular change of variables is no longer required 
(in s > 0); we deduce that the asymptotics of R(^)f have the form 



exp 

(6.9) 



((fci/2 + sk 2 /2 + iy/fa-k*/4)(l + s*j) log a*) (-1/ log:n) 1/2 5 
J l/2 4 2/2 (exp VOoga^ + floga:^) { _ lf l ogXl )V 2 g. 



At the points of the front face of X where log a; 2 / log x% = s this coefficient function 
g is a nonvanishing multiple of 

,s(l + s 2 )-^ 4 (p - fc 2 /4) 3 / 4 5 ( M ?( S )) =a(l + s 2 )- 5/4 ((/i - fc 2 /4)- 3/2 )" 1/2 



xa -*i/a-V^W-fc?/* a .-fa»/a-V/'-MSW-^/4 i2i0i o (a))J2a ^ _ M o (s))/) 
or, finally, 

(6.10) CS (1 + S 2 )- 5 / 4 (/i - fc 2 /4) 3 / 4 F 1 t ( / i?( S ))P*( M - M ?( S ))/. 

Note that although this computation is done separately for each value of s, this 
final expression depends smoothly on s and in fact g is smooth on X up to the 
front face. 

To be precise, we still need to examine what happens near and at the corner 
s = 0, — l/logX2 = 0. However, this is hardly different from the discussion at 
dM\ H Mi. Indeed, we simply introduce t = y/ fj, — Jbi — k 2 /4: as the new smooth 
variable of integration, and then stationary phase can be performed uniformly in s 
down to s = 0, with the limiting contour being the same as for the discussion at 
dM\ n M2, showing that i?(/i)/ is actually poly homogeneous. In particular, note 



that the expansions (6.3), (6.6) and ( |6.10| ) at the front and side faces of X match 
up at the corners, in agreement with the polyhomogeneity. The one point to note 
is that, for example, there is a factor of (— 1/ logxi) 1 / 2 in ( 6.1C ), whereas ((^6) has 



a factor of (—1/ logxi) 3 / 2 . To reconcile this, observe that (3.6) has an extra factor 
of s, and 

K-l/logxi) 1 / 2 = Zj2i£a(_i/i 0ga;i )i/a = (-logx 2 )(-l/logxi) 3 / 2 , 
-logxi 

so in fact the powers match up. 

We remark that if the contour 7 does not go through the critical point /tx?(s), 
then it must necessarily contain points where the integrand is larger; however, at 
those points the phase function itself, not just its imaginary part, will fail to be 
stationary, and so stationary phase gives a decay rate 0((— l/logx)°°); this merely 
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indicates that the contour has not been chosen optimally and a better result is 
possible. 

Before proceeding to the rather more subtle discussion of asymptotics at the 
corner in the presence of bound states, we summarize our results in their absence. 

Proposition 6.1. Suppose f £ C°°(X), p G C\ [fc 2 /4, +00), and Hi, H 2 have no 
L 2 eigenf unctions. Then, with p L = — l/logx;, p _1 = \J p^ 2 + p% 2 , R(p)f has the 
following asymptotic expansion on X: 

(6.11) R(p)f = x\ ll2 x^ 12 eM~Wv~k 2 /4/p)h. 

Here h polyhomogeneous on X, with order 1/2 on the front face and 3/2 on the 
side faces of X . Moreover, the principal symbol of h is an elliptic f -independent 
multiple of 

PM{s))H{H - on the front face, 

(6.12) S 1 {k 2 1 /A)Pl{p - fcf/4)/ on the lift o}M 1 x d~M 2 , 
P{(p - kl/A)S 2 {kt/A)f on the lift of d~Mi x M 2 , 

with Pi(s) given by (|6.§|). 



This discussion already indicates that s plays a dual role, both as a coordinate 
on the space X and also as a spectral parameter. The function Pi(s) identifies the 
front face of X with the line joining the two threshold values k 2 /A and p — k 2 /A 
in the spectral plane. This role becomes even more pronounced in determining 
asymptotics at the front face in the presence of bound states, for then the optimal 
locus for the contour 7 determines which of the residue terms corresponding to the 
different eigenvalues must be included in the asymptotics. In addition, for certain 
real values of p, there are a finite number of exceptional values of s at which p\(s) 
equals either Ah or p — \ 2 i', then the optimal contour must pass through this pole, 
and this creates an additional Legendrian singularity. We explain all of this now. 

6.3. Asymptotics at dM\ x dM 2 in the presence of L 2 eigenvalues for 

p £ C \ M. Thus suppose that Hi or H 2 have L 2 eigenvalues. As we have already 
seen, these correspond to poles of the integrand and the residues at these poles 
may contribute to the asymptotics. To see when this happens, first note that the 
additional terms coming from these residues have the form 



(6.13) 4 2/2+V ^ Al -~ fci/ Vu ® 9u, hi G C°°(M 2 ), 

with analogous terms when the roles of M\ and M 2 have been switched. Here 



The first main observation is that the terms (6.13) dominate the previous one (p.9| 



coming from the continuous spectrum only when 

Im ^Xu-kf/A + sy/ 'p-Xu-k 2 /^ > Im ( ^Jp - k 2 /^l + s 2 ) , 
or in other words when 

(6.14) lmF(X u ) >ImF( M ?(s)). 
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In all other cases, these terms coming from the bound states are lower order in the 
expansion ( |6.9[ ). 

To analyze this, we must examine the geometry of the function F a bit more 
closely Set 

A( M ) - C \ ([fc 2 /4, oo) U (/x - [fcf/4, oo)) , 

and let L C A(/i) be the open line segment connecting the two threshold values 
fc 2 /4 and zx — fef/4; this segment is the image of the map (0, oo) 9s-> Mi( s )- The 
function F depends implicitly on both fj, and s; we shall hold ix fixed throughout, 
but shall study how F changes as s varies. We have already noted that ImF has 
a saddle point at ^i(s). To see this, observe first that this is the unique critical 
point. Next, 



lim lxsxF(fj,i) — s Im \/ fi — fc 2 /4, lim ImF(zxi) = Im \J [i — fc 2 /4, 



and both of these values are greater than \m.F(n\(s)) = \/ 1 + s 2 Im — /c 2 /4 
since Im \J [i — k 2 /4 < 0. Also, each of the square roots in the expression for F has 
imaginary part vanishing along one of the slits and tending to —oo as (i± moves 
horizontally away from that slit, and so ImF tends to — oo as — > oo. Thus 
Hi(s) is indeed a saddle point for ImF. We make the following definition. 

Definition 6.2. For each s, let M s and V s denote the regions in A(/i) where 
ImF(fii) < ImF(^(s)) and ImF(/ii) > ImF(/x°(s)), respectively. 

The union of M s and the slits always has one component, while V s has two 
components, each of which has compact closure intersecting precisely one of the 
two slit axes and which intersect at Hi(s). Label these two components V l s and "PJ, 
respectively, the I and r denoting whether these regions touch the slits extending 
toward the left or right. Finally, recall that the 'optimal' contour of integration 7 
is any contour which remains entirely within the region M s , except when it passes 
through the saddle point. 

Now consider the locations of the poles Ah and \2i relative to the contour 7 and 
these regions M s and V s ■ For any s, these two collections of poles are separated by 



the initial choice of contour 7 in (2.1). However, in deforming 7 to a 7' which lies 
entirely within Af s it may be necessary to shift past some of these poles. This is 
not necessary if all of the Ah lie within J\f s U V r s and all of the fi — \n lie within 
M s U V s . When s is sufficiently small, then in fact all of the Ah lie within M s U "PJ; 
all of the n — \2i lie in this region too, and none lie in V s . On the other hand, 
when s — > 00, both the Ah and the X2i lie within M s U "Pf . As s increases from 
to 00, the region V l s gradually engulfs each of the points Ah, and so the optimal 
contour (which must always lie in M s ) must cross these before this happens. A 
similar phenomenon occurs for the \2i- 

6.4. Asymptotics at dMi x dM2 in the presence of 1? eigenvalues for 

fj, £ K. On the other hand, when \l € R and there are bound states in the interval 
(/1 - fc|/4, fc 2 /4), i.e. when 

fi < min {fc 2 /4 — maxjAi^}, &|/4 — max{A2i}}, 

then this analysis must fail for certain values of s, namely those values where either 
Hi(s) = Ah or /i — /^(s) = \2i' for some . For then an optimal contour would 
need to pass directly through a pole. It is clear that there must be some sort of jump 
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Figure 4. The regions V*, V r s and J\f a for small s > 0. The points 
A, B and C denote [i — A23, M — X22 and — A21 respectively. The 
contour 7 had to be shifted through point C, but not yet past A. 
The point B can made to lie on either side of 7; the residue there 
gives a smaller asymptotic term than the critical point Mi( s ) of F. 



in the behaviour at these points, because for nearby values of s this phenomenon 
does not occur and the previous analysis may be used to get the asymptotics; the 
residue at the pole would be included in these asymptotics if s varies slightly to 
one side, but are not included if it varies slightly to the other. We now prove that 
-R(m) has a Legendrian singularity at this front face of X at those values of s where 
this occurs. 

To be definite, suppose that sq satisfies Mi( s o) — Ah. Let 7 be the original 
contour defining -R(m), as in Figure |l|, and let 7 S be an optimal contour for any 
value s near s - We may arrange that 7 S (0) = Mi( s ), so 7s (0) = An. When s > so 
(closer to the lift of Mi x 9M2) the contour has crossed past Ah, and the residue 
term is present, whereas when s < so (closer to the lift of dM\ x M2) this residue 
term is not included. 

We proceed as follows. First consider what happens when s — > so from the left, 
i.e. with s < sq. Equation ( |6.9[ ) describes the asymptotics when s < Sq, and we are 
interested in the the behavior in the limit. Write 

i? 2 (At ^ Mi) = #2(m - Ah) + (m - An)i? 2 (M - Mi), 
i.e. expand R2 to first order in Taylor series in \i\ around Hi(sq) = Ah. Thus, 

6/ \ #2(M-Mi) --R2O- Ah) 
rl 2 {n - Mi) = \ 

Mi ~ A li 

1 

R' 2 ((r(p - hi) + (1 - ct)(m - Ah)) da, 

and in particular, lim^^x^ i? 2 (M ~ Mi) — "^(m — Ah). 

Our first remark is that the contour in J i?i(/ii)i?2(M — ^u)f ^Mi can be shifted 
farther from the spectrum of Hi since now the R2 term is independent of /ii, so its 
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contribution is negligible. Hence we only need to analyze J Ri{^i){pi—\n)R 2 {iJL— 
fii)dfii. But 

Ri(pi) -ffeG" - Mi) = #i(mi)(mi - Ai,)-R 2 (m - Mi) 

is holomorphic for mi near Ah, so we may move the contour to go through Ah when 
s = sq. Moreover, from the usual parametrix identity ( |3.3| ), using the holomorphy 
of = (mi — Mi)Ri(pi) as a bounded operator on L 2 near mi = Ah, (mi — 

Xii)Ri(m) has the standard asymptotics. Note that i?i(An) = IIh. 
Thus, the holomorphic function -Ri(mi)^ 2 (M2)/ has asymptotics 



5 , n d / w fci/2 fc 2 /2 iy/fJ.i-kl/4 iy/ ^-k^/i . . 

Ri{fJ-i)R2(H2)J = x 1 ' x 2 ' x x y x 2 v a(Mi,M2), 

where a(/ii, M2) is C°° on J7 Ml x U^-^ 2 x Mi x M 2 , and holomorphic in mi and /i2, 
where U is a neighborhood of Ah. Hence 

5 t \5 t \ f Ri(fi,i)Ra(fJ,- m)f - Ri(ni)Ih(fJ>- Ah)/ 

^HMiJ^iM - Mi)/ = t 

Mi — Mi 

fei/2 fei/2 / > s_l / iy/^-kj/4 V^ _ ' I1_fe i/ 4 < \ 

= x 1 'x 2 ' -(Mi -Ah) -l^i z 2 v a(Mi,M-Mi) 



-^1 £ 2 a(Mi,M-Aii) 



ii ii k?/4 i\j [i — [11 — h\ 14 , \ \ • 1 r 

Now add and subtract x x v a; 2 a (Mi; M — Ah) m the numerator 01 

this fraction on the right, and note that the term involving the holomorphic function 

q(Mi,M - Mi) - a(Mi;M - Mi) 
Mi — Ah 

has the same asymptotics as if Ah were not a pole; thus we only need to consider 



7« 



iy/ni-kl/i iy/fJ,-fJ.i-k%/4 iy/fit-kf/4 i^n-\n-k^/4 
X l X 2 ~ X l X 2 1 \ \ 7 

a(Mi,M - Ah) dMi- 

Mi — Mi 

We can also write o(mi, M ~ Ah) = q,(Ah, m — Ah) + (mi — An)ai(Mi, M — Aij), where 
ai is holomorphic in mi and smooth on M 1 x M 2 ; the term involving a\ again yields 
an expression with the standard asymptotics, as if Ah were not a pole. Hence it 
remains to consider 



iy/ni-kj/4 i^J [i-^i-kl/A i^/[ii~k\/A iy/ii-\ li -k%/4 
x l 

This depends on / only via 



I X-i Xo X-i x 2 , 

a(Aii,M-Aii) / dMi 

7s Mi — Mi 



a(An,M - Ah) = -x x v x 2 (U u <g) R 2 (fj, - Mi)) J 

= -»i (IIh ® P|(/x - Ah))/. 

The factor jej fe i/ 4 here corresponds to the precise decay of the eigenfunction 

4>u. The remaining integral is an explicit function which depends only on s, x\ and 
x 2 (and of course Note that in the region x 2 > 0, it is simply a smooth function 
of s, and x 2 (or equivalently, s and x\, keeping in mind that s = log x 2 / log x\ 
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is near so € (0, +00). Thus we only need to understand the asymptotics of this 
function as 2:1,2:2 — > 0. 

Factor out xf^ = x ?0ZM) x WM-fW'))) to get 

(6.16) 



Mi — Ai, 



7s 



where we wrote the integrand in terms of p\ = — 1/ log x± > to keep the signs 
easier to follow. Note that Im(F(pi) — F(p1(s))) < along 7,. 

So far any contour 7 S that stays in Af s has been suitable for our calculations. Now 
we impose an additional condition, namely that there exists a fixed interval [— e, e] 
such that for so — 6 < s < so j s stays in the region where Im y/p — pi — > 
Im y/JT— Aii — fcf/4. Since this region is one side of a parabola, see Figure ||, this 
condition can be easily satisfied. This condition ensures that the real part of the 



exponent in the numerator in (6.16) is non-positive, hence bounded even as pi — ► 0. 



In addition, the fraction in ( |6.16 ) is bounded by Cp 1 1 for some C > 0, since dividing 



the denominator by pi yields a bounded function. 
The key point is that 

F{p x ) F(n1(s)) = ( Ml - pl(s)?F 2 {pi /i°( S ), M °( S )), 

with F2 non-zero at (0, /^(s)), since p\{s) is a non-degenerate critical point of F, 
while 



\J~H - A H - kl/A - xj p - /ii - fcf/4 = (/zi - Ah)-Fi(a*i - Ah) 

with fi nonvanishing and holomorphic near 0. This implies that the first exponent 
is a smooth function of (pi — p\{s)) 2 / pi, with a negative real part of the same order 
of magnitude as p\ — ► 0, hence Schwartz as (/^i — p\{s)) 2 / p\ — > 00, and the second 
one is a smooth function of (pi — X\i)/pi with a negative real part. However, we 
can only make the Schwartz conclusion here if the real part of the same order of 
magnitude as (pi ~ p\- Since 7 S is on one side of a parabola, this is impossible 
to accomplish if j s is to be smooth (for then it is tangent to the parabola). Hence, 
we need to break up "f s into two integrals which will yield boundary terms. These 
boundary terms cancel when s ^ So- 

The geometry is as follows. Start with the space Zo = [0, l) Pl x (so — S, so] s x K{. 
We first blow up t — 0, pi — parabolically (recall that t is the parameter along 7 S ; 
7s(0) = Pi(s)), so that in the interior of the front face s, tj p\j 2 and p\ become valid 
coordinates. Then (F(pi) — F(p®(s)))/pi is smooth near the interior of the front 
face, and e i ' F< - AI1 ' _i? ' All< - s '- ) - ) / pi is smooth on the blown up space with infinite order 
vanishing off the front face. Away from s = s , e _ls( V M - Ml -fc|/4)/ Pl . g 
rapidly decreasing, hence the standard push-forward theorem for polyhomogeneous 
functions [0 yields the same result as stationary phase. In fact, we should think 
of integrating densities, hence we need to change dpi, i.e. dt, to p 1 / 2 d{t/ p 1 ^ 2 ), 
giving the stationary phase asymptotics p 1 ! 2 times a smooth function of (s, p 1 ^ 2 ) 



(multiplied by the exponential x 1 ^^ 1 ^ that we took outside the integral in ( 6.16| )), 
though the terms of the form p times a smooth function of p and s in fact cancel 
(these are the boundary terms of the integration). This is an alternative way 
of thinking about complex stationary phase when the real part of the exponent 
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is non-positive and at least comparable to or larger than the imaginary part, as 
suggested to us by Richard Melrose. It simplifies our previous arguments here, 
although unfortunately it does not directly apply when considering asymptotics at 
the continuous spectrum, as we do in the next section. 

At s — so we also need to resolve the geometry of the second factor in the 



integrand of (3.16). Had we not performed the parabolic blow-up above, it would 
suffice to blow up 

F{ = {(pi, s, t): t = 0, pi = 0, s = s } C Z Q 

in the usual spherical (homogeneous) sense. Indeed, the exponent in the second 
factor as well as the denominator are smooth functions of (pi— \u) / pi, hence of t/pi 
and (s — so)/pi, and the exponential is rapidly decreasing as \t/pi\ + \(s — So)/pi| —> 
oo since 

Imiy/fi-Xu-kl/A- yJp-pi-kl/A) < -C'\p! - A M | < ~C(\t\ + \a- s \). 

Thus, the second factor is polyhomogeneous on this blown up space, of order — 1 
on the front face, order off the front face. 

Unfortunately these two blow ups, i.e. the parabolic one of t = 0, pi = 0, and 
the spherical one of t — 0, p\ — 0, s = so, conflict with each other, and we need 
to find a common resolution. One common resolution is the following: we blow up 

1/2 

the boundary p\ = of Zq by introducing p = p^ as our new boundary defining 
function. We write (^0)1/2 f° r the blown up space. The spherical blow-ups of 

Fi = {(p,s,t) : p = 0, t = 0} and F 2 = {(p,s,t) : p = 0, t = 0, s = s } 

commute with each other since F 2 is a p-submanifold of F%, so [(^0)1/2; Fi, F 2 ] — 
[(^0)1/2; F2, Fi]. Let ffj denote the front face of the blow-up of Fj, j — 1,2. Near 
the interior of the front face of the second blow-up, i.e. ff2, t/p, (s — sq)/p and p 
are valid coordinates. Now blow up the submanifold 

^3 = {p, (s - s Q )/p, t/p) : t/p = 0, (s - so)/ P = 0, p = 0} 

(which is a single point) to obtain 

z = [(^0)1/2! Fi,F 2 ,F 3 ]. 

This introduces, in particular, the front face of the spherical blow up of 

Fi = {(pi,t, s): t = 0, pi - 0, s = s } C Z Q 

in Zq, except that the boundary defining functions differ. In other words, if we 

1 /2 

blow up the front face of [Zq;^] to admit p^ as a smooth function, then a 
neighborhood of the interior of the front face is naturally diffeomorphic to a neigh- 
borhood of the interior of the front face ff 3 of the blow-up of F3. This can be 
seen explicitly since local coordinates which are valid in this region are given by 
(t/p)/p — t/p, ((s — so)/ p)/ p = (s — so)/ p and p = p 1 / 2 . Since upon blowing 
up F 2l the lifts of Fi and F3 are disjoint, their blow-ups commute. Thus, we can 
rewrite Z as [(^0)1/2; F 2 , F3, Fi\. Since the first factor of the integrand of ( |6.16 ) is 



polyhomogeneous on [(Zo) 1 / 2 ; Fi], while the second factor is polyhomogeneous on 
[(Z )i/ 2 ; F 2 , F 3 ], we deduce that the integrand is (one-step) polyhomogeneous on 
Z. 

While the geometry of the integrand only requires these blow-ups, we need fur- 
ther blow-ups to create a b-fibration for the push-forward. Unfortunately, with the 
blow-ups discussed above, all we can hope for is to create a b-fibration with base 



2,x 
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given by a double blow-up of (lo)i/2 — [0, l)p x (so — S,sq] (which is already a 
blow-up, in the sense of change of the C°° structure, of Yq — [0, l) Pl x (s — 5, sq}). 
Namely, one first blows up 

G 2 = {(p,s) : s = s , p = 0}, 

in (io)i/2j an d then 

G 3 = {(p,s) : (s-s )/p = 0, p = 0} 
in [(Yb)i/2j G2]; we denote this space by 
(6.17) Y=[(Y ) 1/2 ;G 2 ,G 3 }. 

It is then straight-forward to carry out appropriate blow-ups in our resolved space 
Z to obtain a new space Z' such that Z' — > Y is a b-fibration. Hence we can 



apply the push- forward theorem of |17J], and deduce that the integral of ( |6 . 16|) 
yields a polyhomogeneous function on Y. More specifically, the result is a one-step 
polyhomogeneous function on Y (keep in mind that p is the boundary defining 
function!) with order 1 on the lift of p = 0, and order on each of the two new 
front faces. 

The asymptotics from s > sq, s — > sq, can be seen similarly. The optimal 
contours in s > sq have been shifted through An, so there is automatically a 
contribution from the pole. We are interested in letting s decrease to So- Since in 

this region the real part of the exponent of x ^v / ' 11_fe i/ 4 m expansion of Ri(pi) 

is less than that of Xi^ 1, , we now need to expand an operator associated to 
Hi rather than one associated to if 2 around Ah. This is a little more delicate, as 
we discuss below. The parametrix identity gives 

Ri{pi) = Pi{pi) - E{pi)Pi{pi) + E{pi)Ri{pi)F{pi), 

with all terms meromorphic. The first two terms are actually holomorphic, so we 
can shift the integral to the optimal location, and even let it go through Ah. Thus, 
we only need to consider the last term. Here the kernel of E(pi) is polyhomogeneous 
on the standard double space Mi x M\, with rapid vanishing on the right face. We 
expand the holomorphic function Ei(fii) as 

E(fn) = Ei(X u ) + (pi - Ah)£i(/xi) 

similarly to the expansion for R 2 {p — pi) above. The contour of the integral 
/ Ei(Xu)Ri(pi)F(pi)R2(p— pi) f dpi can be shifted farther from spec(-Hi), hence 
it is negligible compared to the other terms. So it remains to deal with 

Ei{\u){pi - Xii)Ri(pi)F(pi)R 2 (p - pi)f dpi. 

Is 

But the integrand above again has the standard asymptotics by virtue of the holo- 
morphy of (pi — An)i?i(/Zi), hence completely analogous calculation apply as for 
s < s above. 

We have introduced these new front faces in Yq in order to understand the 
asymptotics in terms of polyhomogeneous expansions. We should certainly discuss 
whether this complicated geometry is simply an artifact of our method, or whether it 
is necessary in the sense that polyhomogeneous asymptotics do not hold in a simpler 
space. We certainly need at least one of the front faces. Indeed, the residue term is 
order in s > so, and is missing in s < sq, where the asymptotic expansion is order 
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1 (in terms of p x ). Thus, it is possible to obtain polyhomogeneous asymptotics 
through s — sq only if there is some blow-up, and hence a new boundary face, 
along which the term of order tends to zero on approach to the lift of s < sq. 
In addition, the exponential in the second factor is smooth only after F[ is blown 
up in Zo, and so it appears that the spherical blow-up of s — sq, p\ = in Yq, 
or at least the presence of the front face of the last bl ow-u p in [(io)i/2) G2', G3], is 
necessary. It then remains to see whether the integral ( 6.1 6| ) is polyhomogeneous on 



[Y ; {(pi, s) : pi = 0, s = s }]. If it is, it must be order on the front face and order 
1/2 on the lift of s < sq. Correspondingly it has order at most 1/2 on the blow-up 
of the corner, which is the interior of the front face [(io)i/2! G2], although with a 
different boundary defining function. Nonetheless, we would expect to see decay 
in the asymptotics at the front face of [(^0)1/2! G2] away from G3, and there is no 



such decay. Indeed, the exponential in the second factor of the integrand in (S.16) 
is rapidly decreasing in the inverse image of this region under the projection, so it 
can be disregarded. Thus the dominant term as p\ — > is an integral of the form 
J s e^i/w^i - Aii) -1 d/ii, a = - iF "(t4( s ))/ 2 > °- Shifting the contour to be 
vertical, the integral becomes i L e~°* / pi (it + (/i?(s) — A^)) -1 dt, with imaginary 
part 

* a + (Ai2(*)-A M )) a<tt - 
Thus, with ^i(s) — Xu = p\^S, S < 0, and changing variables T — t/p\^ 2 1 



T 2 + S 2 



dT < 0, 



independent of pi > 0. Hence there is no decay as p\ — * 0. The limit S — > can 
also be analyzed; indeed, it is the distribution (T + i0) _1 paired with e~ aT . The 
limit from s > sq is similar, but now we get (T — i0) _1 paired with e~ aT . The 
difference, 2iri, is exactly the residue corresponding to the pole. 

Note moreover that the principal symbol of R((J-)f in s < so is an elliptic /- 
independent multiple of P*(/x?(s))P2 (p — Pi(s))f, which in turn is the restriction 
of 



x -U^-i^^ RM{s))p ^ _ M o (fi))/ 

to x\ = 0. That is, the top term of the asymptotics of R(p)f, up to elliptic smooth 
factors and after factoring out x^^x^^e -1 ^/ ^-W/i/p ^ - IS 



^V fcl/2 " VM!w " fc?/4 «i(MSW)^(M-MSW)/, 

which behaves as u , Pl ,_. x 1 kl ^ 2 ^Mit 5 ) fcl/ ' 4 n li P : |(/i — p9(s))f as s — > so (up to 
lower order terms). We have thus explicitly matched the coefficients from the two 
sides of the corner. 

Altogether, we have proved the following theorem. 

Theorem 6.3. Suppose f e C°°(X) ; peC\ spec(PT). Let p, = -l/loga*, p^ 1 = 
y P\ 2 + P2 2 , s — Pi/P2- Then R(p)f has the following asymptotic expansion on 



;so 
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X: 
(6.18) 



R(p)f = x\ l/2 x k 2 2/2 exp(-iy/n ~ k 2 Jl/p)h 



JVi . N 2 , 

Ek 2 /2+iy/ fi-Xn-kl/i . k 1 /2+i v / f j,-\ u -k%/i 

X 2 \4>li ® hli)Xli + 2_^ x l {h>2i ® <P2i)X2i, 

i=l i=l 



where the terms have the following properties. (See Definition 6Jc for the notation, 
Figure /or a picture.) 

(i) fc M G C°°(M 2 ), ft* S C°°(Mi) safe/j/ 

(^h®MIm iX 8M 2 = -(niiOP^CM- Ah))/, resp. 

(6.19) 

(/ijm <8> fe)laM lX M 2 - = -(PfC" - A 2l ) <8> n 2i )/. 

(ii) For p not real, Xji is a function of s, Xji S C°°(R + ) is identically 1 or for 
sufficiently large and small s, with the effect that 

(a) If Xu G V r s , then Xu( s ) — 0, and if \n G V s , then X2i(s) — 0- 

(b) J/Aii G Vj, respectively if X 2l G P J, th en Xu( s ) = 1, r <ssp. X2*(s) = 1. 

(c) Xji G jV s; f/ien the first term in (6.18) dominates the one corresponding 
to <j)ji , hence the choice of Xji is irrelevant. 

(iii) If p is real, p < inf spec(iJ), then xu is a C°° function of (s — So)/ p, so given 
by 



(6.20) 



„ p - k 2 2 /A + slkl/A 
Ah = Mi( s o) = r- — 2 ' 



x' H is compactly supported, limg^+oo Xii (5) = 1, limg^-oo Xii(^) = 0- 5im- 
ilarly, X2% is a function of (s — so)/p, So given by 

(6.21) A 2 i=/x-Mi(so) = r— i > 

1 ' s o 

X 2l «s compactly supported, lim s ^ +00 X2i(5) = 0, limg_ > _ 00 Xh(S) = 1. 

(iv) If p ^ R, then h is polyhomogeneous on X, of order 1/2 on the front face, 
3/2 on the side faces. The principal symbol of h is an elliptic f -independent 
multiple of 

Pi(pi(s))Pl(p: - Pi{s))f on the front face, 

(6.22) Si(fe?/4)P|(> - fef/4)/ on tfie ft/t ofM 1 x <9M 2 , 
PfC/x - kl/A)S 2 {kl/A)f on the lift of dM x x M 2 . 

(v) If p £ R, /i ^ spec(_ff) 7 £/ien g is polyhomogeneous on X doubly blown up at 
the finite number of submanifolds s = sq, p\ = 0, of the front face of X , as 
in ( |6.17|) ; with so given by either of the equations j6.2C| )-( [3.2l| ). The order 
at the front faces of the blow-ups is 0, while the order on the old front face 
is still 1/2. The principal symbol on the front faces depends on f only via 
{Yin® P^p — \u)) f , resp. (Pf (p— A 2 i)<g)II 2 i)/; it is indeed an elliptic multiple 
of these. 
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7. Asymptotics inside the continuous spectrum 

The arguments needed to analyze the asymptotics of R(p ± iO)f when / G 
C° C (X) and p > k 2 /A are not substantially different, and in some senses even 
simpler because the contributions from the L 2 eigenf unctions below the continuous 
spectrum are always negligible away from the side-faces. Since threshold eigenvalues 
introduce additional complications similar to the discussion of the previous section, 
with the additional issue that the push-forward results of JItJ are not directly 
applicable since the stationary phase arguments have a substantially different flavor 
when the phase is pure imaginary, we will assume that k 2 /A is not a pole of Rj. 
Thus, we have the following theorem. 

Theorem 7.1. Suppose f G C°°{X), p > k 2 /A = (k 2 + k 2 )/A, and k 2 /A is not a 



pole of Rj, j = 1,2. Then, with pi = —1/ \ogXi, p 1 = W p x 2 + p 2 2 ' , R(p — iO)f 
has the following asymptotic expansion on X : 

R(p - iO)f = x\ l/2 x k 2 2 ' 2 exp(-i^p-k 2 /A/p)h 

(7 1\ Ni , N 2 , 

+ }^ X 2 (<?li ® hii) + } X x v (h 2 i ® 4>2i), 



i=l 



where h is one-step poly homogeneous on X, and hn and h 2 i are one-step polyho- 
mogeneous on M 2 and M\ respectively. Moreover, h has order 1/2 on the front face 
and 3/2 on the side faces, while the hu is in C° C, {M 2 ) (i.e. has order at dM 2 ) 



and h 2 i £ C°°(Mi). The principal symbol of h is given by (3.22) 



Remark 7. .2. Note that due to the exponential decay of a; 1 (f>u for non-threshold 
eigenvalues \u, namely 



xl^^u G xi kll ^C^{M x \ k 2 /4 - Mi > 0, 

and the similar decay of x 2 k2 ^ 2 (p2i, the second term is dominated by the first one 
everywhere but at the lift of M\ x dM 2 , and similarly the third term is dominated by 
the first one everywhere but at the lift of dM\ x M 2 . In particular, on a neigborhood 
U of the lift of M\ x dM 2 , one has the asymptotics 

R(M)f\v = 4 l/2 4 2/2 expH ^fu~WAlp)g + £ x^^^^^u ® 9u), 

i=i 

while on a neighborhood U' of the front face which is disjoint from dM\ x M 2 and 
Mi x dM 2 both the second and third terms are irrelevant: 

R(j*)f\u> = x kl/2 x^ /2 exp(~i^p-k 2 /4/p)g. 

Remark 7.3. The asymptotics is unchanged away from the corners of X even if k 2 /A 
is a pole of Rj, j = 1,2, in the sense that these do not contribute to the asymptotics 
in the interior of the front face of X, and give the standard contribution in each 
appropriate side face. This follows from stationary phase arguments using the 



integral representation ( 5.15 ). In fac t, we can always reduce the calculation to the 



evaluation of an integral like ( |6.16 ), though here it is conven ient to rewrite the 



fraction appearing in that formula as an integral analogous to ( 6.15 ) 
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Figure 5. Contour of integration 7' to describe R(fi) inside the 
main branch of the continuous spectrum. The cuts C\ and C 2 , 
corresponding to the spectra of H\ and -ff 2 , respectively, have been 
rotated to the cuts C\ and C 2 , as in Figure ^, to make the picture 
clearer, although the analytic continuation is not needed for our 
proof. 

Proof. This is simply a limiting case of the argument of the previous section, but 
which has additional simplifying features. Thus, we simply consider a smooth 
contour 7' that runs along the real axis on the interval fi — fcf/4], but avoids 

the poles of Ri as well as \i minus the poles of R2, see Figure |j| Such a contour 
is simply the limit of the contours 7' that we have considered before. Note, in 
particular, that with 

= yJfH-kl/A + SyJn-m-kZ/4, S = ^i 

as in the previous section, ImF^/ii) < everywhere, and it is strictly negative when 
Hi [fcf/4, /x— fcf/4]. Since the critical point of F lies in [kf/4, — k^/i], the choice 
of the contour outside this interval is irrelevant (except that it should be far from 
the real axis near infinity, as before, to ensure convergence). In particular, Xji G J\f s , 
for all s G [0, +00), except if s = 0, when ImfY/Lt — A 2i ) = (but luiF(Xu) < 0); 
of course a similar statement holds at = with A 2 i replaced by Ah. Thus, 
the residues of the poles of Ri and i? 2 also provide lower order contributions than 
the stationary phase term in most regions in X, except that the poles of Ri give 
the same order as the stationary phase term at Mi x 9Af 2 , and similarly for i? 2 . 
The standard stationary phase lemma now gives the desired results, much as in 
the previous section. To deal with the points \i\ = kf/A and /ii = /J, — fc|/4, i.e. 
the endpoints of [kf/4, /J, — fcf/4], one introduces t% = y//i — kf/4, resp. r 2 = 
yfJL — jLti — fc|/4, just as in the previous section, and notes that in < s < C, 
the phase with respect to t\ is never stationary, while with respect to r 2 it is only 
stationary at s = 0, hence the left end point gives a rapidly decreasing contribution 
in this region, while the right end point gives asymptotically non-trivial terms only 
at s = 0, exactly as expected. □ 

8. The resolvent compactification 

We now turn to a problem originally posited as one of our main goals, namely 
to define a resolution Xj? es of X 2 , which we call the resolvent double space, which 
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is a manifold with corners and on which R(p) is polyhomogeneous, at least when 
p £ spec(-ff). This space is meant to be an analogue of the 0-double space (Mj)f t , 
for either of the conformally compact factors Mj, which we discussed in §3. The 
existence of this space sets the stage for all further development of the analytic 
properties of the Laplacian, including the scattering theory, on X. 

In fact, the last section contains essentially all of the requisite analysis, and what 
remains here is to show how those calculations may be interpreted. As usual, the 



starting point is the integral representation (2.1) for R(p). Again we begin by 
assuming that neither operator Hj has bound states. 



First, as explained in Theorem 3.1, apart from the usual diagonal singularities, 
the Schwartz kernels of both Ri(p\) and R 2 (p — Pi) are polyhomogeneous at the 
boundary hypersurfaces of the 0-double spaces (Mj)q. Hence the Schwartz kernel 
of the integrand i?i(/ii)i? 2 (/i ~ Pi) is polyhomogeneous on the product of these 
0-double spaces, (Mi)q x (M 2 )o, again with diagonal singularities on each factor, 
and thus it is clear that our resolution process should start here. Our task is to see 
what new singularities are introduced by the contour integration, and how these 
may be accomodated geometrically. We shall systematically neglect the diagonal 
singularities in the ensuing discussion; that they produce the correct singularity 
upon integration is obvious in the interior, and requires only minor justification, 
which we omit. 

This product of double spaces has six boundary hypersurfaces, namely 

ffi x(M 2 )o, lfi x(M 2 )l rfx x(M 2 )l (M x )g x ff 2 , (M x )g x lf 2 , (M x )l x rf 2 , 

where ffj is the front face, and If,- and rfj the left and right faces of (Mj)q. The 
kernel of Rjiy) is the product of a smooth function on (Mj)g (albeit with a diagonal 
singularity) with the factor 



0%PrfJ " v 31 ; 

in particular, it is smooth up the ffj. It is straightforward to see that R(p) is 
smooth at the two 'front faces' of the product, fir x(M 2 )g and (Mi)q x ff 2 . This is 
because the exponents in the expansions of R\{p\)R 2 {p — pi) at these faces do not 
depend on p\. It remains to analyze the behaviour at the remaining 'side faces'. 

In fact, the computations of the last section give the asymptotics of R{p) at 
(If 1 x(M 2 )g)U((Mi)gxlf 2 ), and in particular at the corner (If 1 x(M 2 )o)n((Mi)gx 
lf 2 ). To proceed further, the salient observation is that while those stationary 
phase computations were motivated by the geometric picture of letting the variable 
(z\, z 2 ) S Mi x M 2 converge to the boundary (in one of three ways), they can also 
be interpreted purely analytically as a derivation of asymptotics when integrating 
a function of the form x^x^F, where a,\ and a 2 are the appropriate exponents 
depending on p and p\ and F depends smoothly on these variables. What we are 
asking now is the asymptotics of an integral of almost exactly the same form, but 
where each Xj is replaced by pit . p r f . . 

Replace each of the exponential type defining function p at the side faces of these 
0-double spaces by the polynomial type defining function —1/ log/?. In other words, 
logarithmically blow up each of the two side faces If, and rfj of (M j)q, but not the 
front face. Call the resulting space (M,-)q j . For notational convenience, set 



R]i j = ~ 1 / log pif j , Rri j = - 1/ log p Ti j , j = 1 , 2. 
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Now perform the same calculations as in the last section. It follows that the asymp- 
totics of R(p,) along the interiors of each of the four side faces is polyhomogeneous 
in terms of these new defining functions. 

However, to fully resolve the singularities at the various corners, we must perform 
a further sequence of blow-ups. More specifically, the stationary phase calculation 
produces terms polyhomogeneous in R\f . and i? r f . , and the only other terms which 
arise here are polyhomogeneous in (log p\{ 2 p T { 2 ) / (log pif 1 p T { 1 ) . Thus we need to 
find a resolution so that this projective coordinate is C°° when it is bounded, and 
similarly for its inverse. It is not hard to do this. We can certainly assume that 
P\i j ,Pti j < 1/2. Now observe that 

/gjs log(pif 2 p rf2 ) = fogAfg | logPrf 2 

l0g(pif x /3rf i ) l°g Pit i + Plf i fog PU + fog p lfl ' 

which is equal to 

(8 2) — ^ rfl I ^' fl ^ rfl 

Rlf 2 -Rlfi + ^Zrfi Rii 2 -Rlf i + Rri i 

Both terms in this last expression are positive, and hence their sum is bounded if 
and only if both terms are. Clearly we require a space where all of the quotients 

Rlfi/Rlfj, Rlfi/R-cfj, Rrfi/Rrfj i,J=l,2, 

and their inverses are smooth (when bounded). 
Thus first blow up 

(8-3) (lfrnrfx) x (M 2 )* log , x (lf 2 nrf 2 ). 

This may be done in either order since these submanifolds are transverse. Next, we 
must blow up the collection of submanifolds covering 

(If! x(M 2 )l log ) n ((Mi)^ g x lf 2 ), (lfx x(M 2 )g >log ) n ((M^j og x rf 2 ), 

(rf! x(M 2 )2 Jog ) n ((Mi)g, log x lf 2 ), (rfi x(M 2 )* log ) f] ((M^ x rf 2 ). 

These are disjoint now after the initial blow-up, and so may be blown up in any 
order. The manifold X^ cs obtained in this way has all the desired properties. In 



fact, working in a region where (8.1) is bounded, we must show that both terms in 
( fT^ ) are smooth. But 

-Rifi Riix _ Rif! 1 



i?if 2 R ifl + R lfl Ru 2 (Rif ! /Rrh ) + 1 : 



with a similar expression for the other term. The first spherical blow up ( 3.2 ) makes 
Rrf./Rif. and its inverse smooth where bounded, while the prefactors Rn 1 /R\f 2 
and R\f 1 /R I { 2 also become smooth because (If i x (M 2 )q log ) n ((Mi)q log x lf 2 ) and 
(If i x (M 2 )q i og )n ((Mi)q log x rf 2 ) have been blown up. We thus conclude that when 
Hi and H 2 have no L 2 eigenfunctions, R(p) is polyhomogeneous on this space X 2 CS . 

To conclude this discussion, suppose that cither Hi or H 2 have bound states. 
Then R{p) will have additional singularities at certain submanifolds of X 2 CS where 
log(pif 2 /o r f 2 )/log(pif 1 p r f 1 ) assumes the same cri tical values sq as already appeared 
for s = log x 2 / logs i in part (v) of Theorem 3.2 . That is, these singularities happen 
at the submanifolds given by S — sq inside dX 2 cs , where S is the smooth function 

° n X L g iven b y 

* = T, 1" T> — 7~F> T~F> — = { a lU2 + tT Hr2) 



R\f 2 Rrf 2 J {Rlfi/Riii) + 1 T Carl + l' 
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in the region where aim = "s^i ai\ri = and cziri = RxuIRtU are valid 
coordinates, i.e. where they are bounded. Since the second factor is thus bounded 
from both below and from above, we conclude that, as long as sq — is not one 
of the critical values, at least one of 07102 and aii r 2 is non-zero at any point on 
S = so, and in addition d aill2 S ^ 0, d ailr2 S ^ there, so S — sq is a codimension 
one p-submanifold of X 2 CS , and we have singularities at its intersection with dX 2 csl 
which in turn is a finite union of codimension one p-submanifolds of dX 2 cs . These 
singularities can be resolved by the same double blow up, preceeded by the square 
root blow up of the boundary defining function, as on the single space. 

While these blow ups may appear somewhat complicated, the only relevant part 
of X 2 0S as far as the Martin boundary is concerned, is the inverse image U' of 
regions X x U, U compact subset oi X , U = Ui x U2 C M\ x A/ 2 open, under the 
blow-down map f3 : X 2 CS — > X 2 . Thus, in the second factor we always stay away 
from the boundary, hence any blow-ups involving rfi, rf2 can be neglected. Thus, 
the intersection of U' with X 2 CS is a subset of 

[(Mi)iog x U x x (M 2 )io g x U 2 ;d(M 1 ) log x U x x d(M 2 ) log x U 2 ] 

which is essentially the same as X x U\ x U2 (i.e. they agree up to the rearranging 
of factors). Hence, in this set the asymptotics of the kernel of R(fi) is given by 
Theorem |6.3| , with an extra variable on U added, and references to / dropped (this 
really corresponds to applying the resolvent to delta distributions at points z and 
letting z vary). 



9. The Martin compactification 

With the information we have collected and proved about the resolvent R(fJ.), it is 
now a simple matter to determine the full Martin compactification of X = Mi x M2. 
We refer back to §5 for the general details of how this construction is to be carried 
out, but we briefly recall the main ideas. Fix any /1 G M, fi < infspec(iZ). In 
addition, let p S X be fixed and = {w[ e \ u4 ) € X any sequence which leaves 
every compact set. Define 

H{fi;p, w (l >) 

this function is a solution of (H — n)U e = for z ^= and is normalized so that 
U^\p) = 1. Any such sequence has a subsequence U l which converges uniformly 
on compact sets to a function U which satisfies (H — n)U — on all of X and 
U(p) = 1. The Martin boundary Xm is defined to be the set of all possible limit 
eigenfunctions that arise in this manner. 

First consider the case if neith er H i nor H2 has L 2 eigenvalues. In this case, 
by the asymptotics of Proposition ^]l] (see also the last paragraph of the previous 
section) 

(9.1) RMlxxx = xi' 2 4 2 ' 2 exp(-VM - fc a /4/p)M, 

where ho is the product of the square root of a defining function of the front face 
of X with the 3/2 power of a defining function of each side face of X, and g is C°° 



Mi 
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in X x X with restriction to each face being an elliptic multiple of 
Plfa°i(s))P^fa - n°(s)) on the front face, 

(9.2) Si(k\/A)P^fa - fcf/4) on the lift of Mi x d~M 2 x U, 
Pi fa - kl/A)S 2 (kj/4) on the lift of dMi xM 2 xU. 

In particular, 

Rfa;z,w) 

Ufa,z,w) = — - 

Rfa;p,w) 

is a continuous function on X w x X z . Thus, the map 

(9.3) Ufa, ., io) € C°°{X) 

defines a map X — *• dX m which is continuous and surjective. More explicitly, 
if w^' is any sequence of points in X converging to w G dX, the continuity of U 
implies that U^\z) = Ufa,z,w), hence Ufa, .,w) is a point in the Martin boundary. 
Conversely, suppose that w^' is any sequence of points in X which leaves every 
compact subset of X, and such that converges uniformly on compact subsets 
of X . By passing to a subsequence ' we may assume that w^ e ) converges to 
some w £ dX since X is compact. Hence ftmU^ \z) = Ufa,w,z), so the map is 
indeed surjective. 

Explicitly, U is given by 

(9.4) 

Ptfa\{s)- Zl , Wl )Plfa-^{ S )- Z2 , W2 ) {s ^ w2) e { ^ +ooh x ^ x ) 



Pffa° 1 (s);p 1 ,w 1 )P 2 t fa - fi1(s);p 2 ,w 2 )' 
S 1 (k 2 1 /A;z 1 ,w 1 )P 2 t fa-k 2 2 /A;z 2 ,w 2 ) 
5i(A?/4;pi,ici)i^(/i- kj/4;p 2 ,w 2 )' 
Pi fa - fcj/4; zi, Wl )S 2 {kl/A- z 2 ,w 2 ) 
Pi fa- kl/4;p 1 ,w 1 )S 2 (kj/4:;p 2 ,w 2 )' 



(wi,w 2 ) € Mi x dM 2 , 
(w!,w 2 ) £ dMi x M 2 . 



Thus, the injectivity of the map (9.3) depends on whether these are all different 
elements of C°°(X Z ) as w varies. It is easy to see that the restriction of this map to 
the front face is indeed injective since Pjfa; Zj, Wj) has a pole when Zj — ► Wj and is 
otherwise bounded, and so the points w± and w 2 are uniquely identified, as is the 
eigenparameter fj.i(s) and hence s itself. On the other hand, it is not obviously true 
that the generalized spherical function 5j(fc|/4; Zj, w®) determines the point w®. In 
case Mj is hyperbolic space, this is holds because of the rotational symmetry around 
W® , and so must also hold for conformally compact metrics (globally) near to the 
hyperbolic one. In general, these portions of the Martin boundary are identified 
with Xi x M 2 and Mi x X 2 , where 

Xj = image(w° — > Sj(k^/4; Zj,Wj)). 

Next suppose that Hi has L 2 eigenvalues. The asymptotics of the resolvent 



kernel in X x X are now given by Theorem 3.3, again in the sense discussed in 
the last paragraph of the previous section. Let An be the bottom eigenvalue of Hi 
which is hence simple, and let sq be the corresponding value of s = pi/p 2 given by 
( HI ). In the region s < so, in particular near dMi x M 2 x X , Rfa) has the same 



asymptotics as before, while on the front face of the double blow up of s = so as 
well as in s > so the leading part of the asymptotics is given by elliptic multiples 
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of x 1 kl ^ 2 *v /Ali fel / 4 n n x P 2 {ft — which are independent of z e X. We again 
set 

R(fj,;z,w) 

U(fi,z,w) = — -, 

R{^P, w) 

and U is again continuous onXxX. Now, however, w i— > J7(/i, ., w) factors through 
the map <9X — » <9X C , where dX c is the collapsed boundary of X given by 

dX c = (dMi x M 2 ) U ([0, s ) s x 9Mi x d~M 2 ) U <9M 2 , 

and the collapse map <9X — > <9X C is the identity on the first two sets on the right 
hand side and is the projection [s , +00) x dM\ x dM 2 — ► 9M 2 , resp. Mi x 9M 2 — ► 
<9M 2 in the other parts of dX. Thus, we obtain a continuous surjection from dX c 
to BXm- 

Again, U can be written down explicitly. It is the same as in the case without 
L 2 eigenfunctions when s < s , and it is 

P^(/j,-k^/4;z 2 ,w 2 ) 



w 2 e dM 2 . 



Thus, the map dX c — » 3Xm is injective over the collapsed part of the boundary, as 
can be seen by letting z 2 — > w 2 , and otherwise we are exactly in the same situation 
as in the case where there are no L 2 eigenfunctions. 

Similar arguments work if H 2 has eigenvalues, or even if both H\ and H 2 do. In 
the latter case X collapses near both side faces. 
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